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A  formalism  is  developed  for  studying  the  effects  of  iso- 
topic  impurities  upon  the  conduction  of  energy  in  a  three-dimension- 
al dielectric.   The  model  used  is  the  analogue  of  a  one-dimensional 
model  introduced  by  Allen  and  Ford  in  I968.   The  Kubo  method  of  cor- 
relation functions  is  the  particular  technique  employed  to  determine 
the  current  conducted  through  the  solid.   Emphasis  is  placed  upon 
obtaining  relevant  mathematical  techniques  for  treating  the  problem 
in  three  dimensions. 

Two  special  distributions  of  impurities  within  the  solid 
are  considered  and  solutions  are  found  for  the  current  in  each  case. 
A  detailed  treatment  is  given  for  the  case  of  isotopes  uniformly 
distributed  among  the  host  atoms ,  taking'  into  account  the  effects 
of  multiple  scattering  of  phonons  by  impurities.   A  calculation  in- 
volving the  first  Born  approximation  is  undertaken  for  the  case  in 
which  isotopes  are  distributed  at  random  throughout  the  crystal. 

vi 


The  results  are  compared  with  the  results  obtained  in  one 
dimension.   Some  discussion  of  problems  remaining  to  be  solved  in 
three  dimensions  is  given. 
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CHAPTER  I 
INTRODUCTION 

Attempts  to  place  the  theory  of  energy  transport  in  solid 
insulators  on  a  firm  foundation  and  to  obtain  a  reliable  model  from 
which  calculations  can  be  made  have  been  underway  since  the  early 
1900' s.   Although  research  on  this  problem  has  been  extensive,  the 
inherent  difficulty  in  treating  non-equilibrium  systems  has  posed  a 
formidable  barrier  to  success.   Indeed,  the  very  early  investigators 
were  forced  to  content  themselves  with  qualitative  predictions  or,  at 
best,  calculations  involving  rather  drastic  approximations.   Before 
discussing  the  specific  aims  of  this  dissertation,  it  therefore  seems 
advantageous  to  present  a  historical  survey  of  the  research  done  on 
this  problem  up  until  the  present  time.   The  intent  will  be  to  indi- 
cate how  a  more  quantitative  understanding  of  this  problem  has 
emerged.   It  should  be  emphasized,  however,  that  this  is  not  intended 
as  a  detailed  review,  but  rather  a  brief  and  completely  qualitative 
synopsis  of  what  the  author  considers  some  of  the  more  significant 
advances.   Particular  emphasis  is  given  to  the  work  of  the  last  few 
years  and,  where  applicable,  pertinent  experimental  evidence  has  been 
indicated. 

1.   Phenomenological  and  Related  Theory 

Energy  transport  in  solids  is  most  conveniently  discussed 
by  defining  a  quantity  known  as  the  thermal  conductivity.   From 


experimental  evidence,  one  knows  that  if  a  temperature  gradient 
3T/3z  is  introduced  across  an  isotropic  solid,  the  energy  which 
flows  through  the  solid  obeys  the  Fourier  heat  law  given  by*- 

J   =  -  K  1^    .  I-l 

Z  dZ 

J   is  the  current  or  energy  per  unit  area  per  unit  time  which  flows 
z 

in  the  direction  of  the  gradient  and  K  is  the  thermal  conductivity. 
The  negative  sign  indicates  that  the  current  flows  from  higher  to 

lower  temperatures. 

The  first  evidence  of  any  experimental  effort  to  deter- 
mine the  factors  upon  which  the  thermal  conductivity  depends  is  the 

2 
early  investigations  of  Eucken.    Eucken  found  that  for  most  solid 

dielectrics  the  thermal  conductivity  varied  as  1/T  at  high  tempera- 
tures.   Debye,   later  in  1914,  succeeded  in  explaining  the  1/T  de- 
pendence using  kinetic  theory  and  thermodynamic  considerations.   The 
argument  was  entirely  classical,  not  intended  to  be  rigorous,  and 
provided  no  adequate  means  for  determining  the  proportionality  con- 
stant. 

The  first  really  significant  advance  was  suggested  in  1929 

4 
when  Peierls,   using  the  newly  developed  idea  of  quantized  lattice 

waves,  succeeded  in  formulating  a  theory  valid  at  low  temperatures. 
Peierls'  formulation  was  based  on  the  Boltzmann  transport  equation  and 
formally  included  the  effects  of  the  so-called  phonon-phonon  inter- 
actions.  These  interactions  arose  from  the  retention  of  terms  of 
higher  order  than  quadratic  in  the  crystalline  potential;  they  con- 
sisted of  the  normal  processes  which  conserve  the  crystal  momentum 


and  the  umklapp  processes  which  do  not.   In  effect,  Peierls  was  able 
to  obtain  an  integral  equation  which  represented  a  solution  to  the 
collision  term  in  the  Boltzmann  equation,  though  he  was  unable  to  ob- 
tain a  solution.   However,  his  theory  did  predict  that  the  thermal 
conductivity  should  increase  exponentially  with  decreasing  tempera- 
ture in  the  low  temperature  limit. 

De  Haas  and  Biermasz    attempted  experimental  verifica- 
tion of  Peierls'  results.   They  found  that  the  thermal  conductivity 
did  increase  with  decreasing  temperature  though  not  so  rapidly  as 
Peierls'  theory  had  predicted.   Though  it  was  not  understood  at  the 
time,  the  discrepancy  is  now  known  to  have  arisen  from  the  resistance 
to  the  flow  of  energy  caused  by  point  defects  which  had  not  been  in- 
cluded in  Peierls'  analysis. 

During  the  next  several  years,  a  number  of  noteworthy 
attempts  were  made  to  extend  Peierls'  original  theory.   Specifically, 
investigators  sought  to  account  for  the  various  mechanisms  which 
cause  resistance  to  the  flow  of  energy  in  a  crystal  lattice.   Most 
calculations  involved  a  relaxation  time  approximation  to  the 

Q 

Boltzmann  equation.   In  1938,  Casimir  suggested  that  the  mean  free 
paths  of  phonons  should  be  limited  by  the  finite  size  of  the  sample. 
This  effect  was  known  as  boundary  scattering  and  was  shown  to  be  the 
dominant  form  of  resistance  to  the  flow  of  energy  at  very  low  tempera- 
tures.  Casimir's  theory  predicted  that  at  low  temperatures  the  ther- 

3 
mal  conductivity  should  vary  as  T  and  depend  on  the  size  of  the  spec- 
imen, an  effect  which  had  been  observed  in  the  earlier  experiments 
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of  de  Haas  and  Biermasz.   Later  in  1942,  Pomeranchuk  pointed  out 

that  point  defects,  or  variations  in  the  mass  distribution  through- 


out  the  lattice,  could  cause  thermal  resistance.   As  indicated  pre- 
viously, this  accounted  for  the  fact  that  the  thermal  conductivity 
did  not  rise  exponentially  with  decreasing  temperature  at  low  temp- 
eratures.  Some  years  later,  Klemens   showed  that  the  normal  pro- 
cesses were  unimportant  except  at  low  frequencies  and  that  the  um- 
klapp  processes  were  the  only  anharmonic  processes  which  resisted 
the  flow  of  heat.   Finally,  in  1959,  Callaway,    using  frequency- 
dependent  relaxation  times,  included  the  effects  of  the  normal  pro- 
cesses at  low  frequencies. 

This  technique  of  using  the  Boltzmann  equation  to  inves- 
tigate energy  transport  which  was  begun  by  Peierls  and  extended  by 

the  workers  noted  above  has  come  to  be  known  as  the  phenomenologi- 

-,  ui  12,13 
cal  theory.   Excellent  reviews  of  this  theory  are  available 

and  a  detailed  discussion  will  not  be  given  here.   It  is  interest- 
ing, however,  to  point  out  two  important  consequences  of  the  theory. 
First,  calculations  generally  involve  making  a  priori  assumptions 
concerning  the  manner  in  which  the  system  under  investigation  can 
be  expected  to  return  to  equilibrium  when  disturbed  from  equilib- 
rium.  Second,  the  theory  predicts  an  infinite  thermal  conductivity 
for  those  systems  which  contain  no  scattering  mechanism  for  causing 
their  return  to  equilibrium.   This  second  fact  was  explained  by 
arguing  that  in  any  realistic  crystal  the  divergence  is  eliminated 
by  boundary  scattering  at  low  temperatures  and  by  anharmonicities 

at  high  temperatures. 

14 

Some  excellent  experiments  conducted  by  Berman  and  Brock 

in  1965,  attempted  to  verify  the  results  of  the  phenomenological 


theory.   Fairly  good  agreement  with  theoretical  predictions  was  ob- 
tained for  most  of  the  temperature  ranges  considered. 

We  may  summarize  the  phenomenological  theory  by  noting 
that  the  theory  did  produce  results  that  could  be  fitted  to  the 
experimental  data  reasonably  well.   However,  from  a  theoretical 
viewpoint,  the  theory  was  rather  unsatisfactory.   In  particular, 
the  results  obtained  for  systems  which  contained  no  scattering 
mechanism  were  meaningless.   Clearly  there  was  not  contained  within 
the  theory  a  systematic  method  of  examining  how  the  effects  of 
phonon  scattering  become  important  and  might  be  expected  to  lead  to 
a  theoretical  prediction  of  the  Fourier  heat  law.   Some  more  specif- 
ic faults  of  the  theory  will  be  noted  in  the  following  section. 

In  the  1950' s,  Kubo   developed  a  technique  for  treating 
transport  problems  which  was  not  based  on  the  Boltzmann  equation 
but  rather  which  depended  on  the  solution  of  the  dynamical  equations 

of  motion  of  the  crystal.   Considerable  credence  was  afforded  the 

1  fi 
Kubo  method  by  the  subsequent  work  of  Hardy  and  his  co-workers   who 

managed  to  obtain  from  the  Kubo  formalism  a  perturbation  expansion 
for  the  lattice  thermal  conductivity.   Though  the  analysis  was  rigor- 
ous and  did  succeed  in  placing  the  phenomenological  theory  on  a  more 
satisfactory  basis,  it  was  quite  formal  and  not  readily  amenable  to 
calculations.   In  practice,  it  was  necessary  to  resort  to  approxi- 
mations  which  just  reproduced  the  transport  equation   '    and  which 
led  to  the  same  results  as  the  phenomenological  theory.   Furthermore, 
no  attempt  was  made  to  perform  calculations  using  a  specific  model. 


2.   The  One-Dimensional  Model 

19-21 
During  the  1960's,  several  investigators      independently 

examined  the  possibility  of  using  one-dimensional  chains  of  atoms  to 

study  energy  transport.   Since  the  calculations  of  this  dissertation 

are  based  on  an  extension  of  the  specific  model  used  by  Allen  and 

19 
Ford,    a  brief  discussion  of  the  model  will  be  given.   The  model 

consisted  of  an  infinite  chain  of  atoms  connected  by  Hookeian  springs. 
A  finite  section  of  the  chain  of  length  L_  contained  both  host  atoms 
as  well  as  isotopic  impurities  randomly  distributed  among  the  host 
atoms.   The  infinite  regions  on  either  side  of  this  finite  section 
consisted  only  of  host  atoms  and  were  used  to  simulate  high  and  low 
temperature  reservoirs  which  caused  energy  to  flow  through  the  length 
L„.   A  first  principles  calculation  using  the  Kubo  theory  was  under- 
taken to  calculate  the  current  and  consequently  the  thermal  conductiv- 
ity.  The  model  was  amenable  to  a  formally  exact  solution  in  the  har- 
monic approximation,  but  approximations  were  necessary  to  permit  ex- 
plicit evaluation  of  the  thermal  conductivity.   In  essence,  Allen 
and  Ford  resorted  to  a  Born  approximation  solution  to  the  scattering 
problem. 

It  was  found  that  the  model  produced  a  non-divergent  cur- 
rent even  in  the  harmonic  approximation  and  even  for  the  special  case 
of  an  isotopically  pure  chain  for  which  there  was  no  scattering  mech- 
anism.  Allen  and  Ford  emphasized  that  the  current  should  ultimately 
be  limited  by  the  amount  of  energy  available  for  transport  and  hence 

always  non-divergent.   This  plausible  argument  had  previously  been 

22 
stressed  by  Erdos   who  had  performed  calculations  using  an  indepen- 


dent  method  and  found  a  non-divergent  result  for  the  isotopically 
pure  three-dimensional  crystal. 

Nonetheless,  the  non-divergent  result  was  at  variance  with 
the  phenomenological  theory  and  consequently  Allen  and  Ford  carefully 
reviewed  that  theory.   They  found  that  some  of  the  derivations  in 
the  phenomenological  theory  were  not  internally  consistent.   Specif- 
ically, it  was  noted  that  some  of  the  approximations  led  to  the  pre- 
diction of  negative  occupation  numbers  of  the  phonon  states  at  low 
frequencies,  a  prediction  which  clearly  made  no  physical  sense.   It 
is  now  believed,  therefore,  that  inconsistencies  in  the  phenomeno- 
logical theory  produced  the  low  frequency  divergence  in  those  systems 
which  had  no  tendency  to  equilibrate. 

Though  a  divergent  result  was  not  obtained,  the  approxima- 
tions did  fail  to  produce  a  thermal  conductivity  in  the  usual  sense. 
In  particular,  it  was  found  that  the  heat  current  was  proportional 
to  the  temperature  difference  between  the  reservoirs  rather  than  to 
the  temperature  gradient  across  the  length  L^.   This  is  equivalent 
to  obtaining  a  thermal  conductivity  dependent  upon  the  length  of  the 
sample  which  is  experimentally  known  not  to  be  the  case.   Allen  and 
Ford  speculated  that  either  the  multiple  scattering  effects  (omitted 
due  to  the  Born  approximation)  or  the  anharmonic  processes  might 
cause  sufficient  resistance  to  the  flow  of  energy  to  produce  a  length 

independent  thermal  conductivity. 

23 

Early  in  1972,  Greer  and  Ruben   succeeded  in  calculating 

the  multiple  scattering  effects  in  the  one-dimensional  problem  and 
concluded  that  these  were  insufficient  to  produce  a  length-indepen- 
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dent  result.   Some  recent  work  done  by  Allen,  Ott  and  Powell, 

however,  indicates  that  multiple  scattering  is  sufficient  if  one 
assumes  an  appropriate  distribution  of  impurities  within  the  impu- 
rity bearing  segment  of  the  chain. 

The  primary  conclusion  that  can  be  dra\^m  from  these  cal- 
culations made  in  one  dimension  is  that,  unlike  the  phenomenological 
theory,  a  meaningful  thermal  conductivity  can  be  obtained  even  for 
the  harmonic  system  with  no  mechanism  for  causing  thermal  resistance. 
In  addition,  the  one-dimensional  model  has  shown  considerable  prom- 
ise of  providing  a  systematic  method  whereby  one  can  study  the  ef- 
fect on  the  thermal  conductivity  as  various  scattering  processes  are 
included.   Ultimately  one  hopes  to  understand  the  manner  in  which 
these  scattering  processes  produce  a  current  proportional  to  the 
temperature  gradient  rather  than  to  the  temperature  difference.   In 
the  discussion  which  follows,  it  will  be  seen  that  similar  conclu- 
sions are  possible  for  a  three-dimensional  model. 

3.   The  Present  Calculation 

For  reasons  indicated  in  the  last  section,  Allen  and  Ford 
have  felt  obligated  to  attempt  an  extension  of  the  one-dimensional 
model  to  three  dimensions.   In  1969,  they  succeeded  in  calculating 

the  thermal  conductivity  for  the  three-dimensional,  isotopically  pure 

25 
crystal,    and  concluded  that  it  was  possible  to  obtain  a  meaningful 

result  in  three  dimensions  for  the  case  in  which  no  scattering  mech- 
anism X'/as  present.   It  will  be  the  aim  of  this  dissertation  to  pro- 
vide a  formalism  whereby  one  can  study  how  scattering  becomes  impor- 


tant  in  three  dimensions,  to  introduce  some  techniques  applicable  to 
calculations,  and  finally  to  perform  some  calculations  for  some  rather 
special  cases.   It  should  be  emphasized,  however,  that  the  aim  is  not 
a  complete  solution  of  the  three-dimensional  problem;  unfortunately, 
the  problem  in  three  dimensions  has  been  sufficiently  complicated  that 
much  work  remains  to  be  done. 

The  model  to  be  used  is  the  analogue  of  that  used  by  Allen 
and  Ford  in  one  dimension.   The  crystal  is  taken  to  be  infinite  and 
contain  a  slab  of  infinite  cross-sectional  area  but  finite  thickness 
L_  along  the  z  axis.  (For  orientation  of  the  coordinate  system,  see 
fig.  1.)   Isotopic  impurities  of  mass  M'  as  well  as  host  atoms  of 
mass  M  are  distributed  throughout  the  slab  of  thickness  L.  which  will 
henceforth  be  referred  to  as  the  sample.   The  infinite  isotopically 
pure  regions  on  either  side  of  the  sample  are  used  to  simulate  high 
and  low  temperature  reservoirs.   They  are  maintained  at  different 
temperatures  T^  and  T  where  T^  >  T„.   The  resulting  temperature  dif- 
ference causes  a  current  to  flow  through  the  sample  parallel  to  the 
z  axis.   This  current,  and  finally  the  thermal  conductivity,  will  be 
calculated  using  the  Kubo  formalism. 

The  restrictions  noted  in  the  one-dimensional  model  are 
also  present  in  this  calculation.  Namely,  only  the  harmonic  term  is 
retained  in  the  crystalline  potential  and  therefore  both  umklapp  and 
normal  processes  which  are  anharmonic  have  been  neglec'"ed.  However, 
these  should  make  a  negligible  contribution  to  the  heat  current  at 
the  low  temperatures  which  will  be  considered.  In  specific  calcula- 
tions it  is  further  assumed  that  the  crystal  has  simple  cubic  sym- 
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Figure  1 

The  model.   Isotopic  impurities  are  distributed  through- 
out the  sample  which  is  labeled  II.   Regions  I  and  III  simulate 
high  and  low  temperature  reservoirs  maintained  at  temperatures  T^ 
and  T„.   The  sample  and  the  reservoirs  are  assumed  to  have  infinite 
cross-sectional  area.   For  purposes  of  calculation  the  y  axis  is 
directed  vertically  upv/ard  and  the  x  axis  into  the  page. 
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metry,  that  it  is  isotropic,  and  that  the  acoustic  approximation  is 
valid.   Finally,  it  is  assiimed  that  all  three  polarization  modes  in 
the  crystal  propagate  at  the  same  velocity.   Though  anharmonicities 
have  not  yet  been  included  in  either  the  one-dimensional  or  the  three- 
dimensional  model,  the  additional  assumptions  made  here  are  merely 
for  the  purpose  of  simplifying  the  calculations.   The  model  does  not 
restrict  one  to  make  these  assumptions,  but  the  author  feels  that 
they  do  not  sacrifice  any  essential  understanding  of  the  problem. 

In  Chapter  II  of  this  dissertation,  the  formalism  necessary 
for  solving  the  problem  is  introduced  and  the  formalism  is  cast  into 
a  form  amenable  to  calculations.   In  Chapter  III,  a  solution  to  the 
scattering  problem  is  obtained  for  the  case  of  a  uniform  distribution 
of  impurities.   Multiple  scattering  effects  are  included.   Attention 
is  given  to  the  details  of  the  calculation  since  the  author  believes 
that  the  same  techniques  may  very  v;ell  be  applicable  to  more  general 
distributions  once  the  mathematical  difficulties  have  been  overcome. 
In  Chapter  IV,  the  heat  current  is  calculated  for  the  case  discussed 
in  Chapter  III.   In  addition,  a  Born  approximation  calculation  of  the 
current  is  made  for  a  random  distribution  of  impurities.   Finally, 
in  Chapter  V,  the  results  are  discussed  and  analogies  with  the  one- 
dimensional  problem  are  noted.   The  relevance  of  the  solutions  ob- 
tained to  the  experimentally  observed  phenomenon  of  Kapitza  resist- 
ance is  briefly  mentioned.   Finally,  some  discussion  of  work  remain- 
ing to  be  done  is  given. 


CHAPTER  II 
FORMALISM 

The  purpose  of  this  chapter  will  be  to  present  the  formal- 
ism necessary  for  calculating  the  current  for  the  model  discussed. 
In  particular,  the  Kubo  formula  for  the  heat  current  will  be  given 
and  a  method  for  obtaining  the  current  operator  which  must  be  used 
in  conjunction  with  that  formula  is  outlined.   Third,  the  Green's 
function  technique  for  obtaining  the  equations  of  motion  of  the  crys- 
tal  is  discussed.   Finally  a  formal  expression  for  the  current  in 
terms  of  the  quantities  to  be  calculated  will  be  obtained.   Second- 
arily, this  chapter  v;ill  aim  to  set  the  notation  used  in  the  remain- 
der of  the  discussion  and  to  present  some  rudimentary  ideas  relating 
to  the  phonon  representation.   In  cases  where  certain  aspects  of 
the  formalism  have  been  developed  elsewhere,  only  results  are  given, 
though  specific  references  are  always  indicated. 

1.   Kubo's  Formula 

The  Kubo  theory  of  correlation  functions  is  a  technique 
which  is  capable  of  producing  expressions  for  the  density  matrix  for 

certain  non-equilibrium  systems.   The  mathematical  details  of  the 

9fi  97 
theory  have  been  discussed  elsewhere   '    and  will  not  be  discussed 

here.   The  general  assumptions  of  the  theory  as  discussed  by  Mori 

et  al.   will  be  briefly  enumerated,  however. 
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In  the  Kubo  theory  it  is  assumed  that  the  system  under  in- 
vestigation may  be  divided  into  macroscopically  small  segments  each 
characterized  by  a  local  temperature  T(r).   Initially,  these  seg- 
ments are  not  allowed  to  interact  with  one  another  and  each  is  as- 
sumed to  be  in  thermal  equilibrium.   At  time  t  =  0,  the  segments  are 
allowed  to  interact  and  bring  about  a  change  in  the  density  matrix 
of  the  system.   It  is  further  assumed  that  the  system  will  approach 
a  steady  state  and  that  the  temperature  as  a  function  of  position 
does  not  change  as  the  steady  state  is  approached.   The  implications 
of  this  last  assumption  will  be  further  discussed  in  Chapter  V. 

At  time  t  =  0,  the  density  matrix  of  the  system  is  given 
by  the  local  equilibrium  distribution 


p(0)  =  I"  exp[-/6(?)  H(?)  d?]  II-l 


where  Z  is  the  partition  function,  H(r)  the  Hamiltonian  density,  and 

->■        -> 
3(r)  =  k  T(r)  where  k  is  Boltzmann  s  constant.   The  time  evolvement 

of  the  density  matrix  is  determined  by  the  quantum  mechanical  equi- 
valent to  the  Liouville  equation  given  by 


where 


in\p^=  [H,p(t)]  II-2 


p(t)  =  exp  (-iHt/Ti)  p(0)  exp  (iHt/Ti)   .      II-3 


In  these  equations,  H  is  the  Hamiltonian  of  the  system.   It  is  de- 
sired to  evaluate  the  density  matrix  at  some  time  T  after  the  steady 
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state  has  been  reached  but  before  the  system  approaches  complete 
thermal  equilibrium.   Thus  for  the  model  at  hand,  we  wish  to  choose 
T  large  enough  so  that  ample  time  has  been  allowed  for  energy  to 
propagate  through  the  sample,  but  small  enough  so  that  energy  has 
not  propagated  all  the  way  out  to  the  ends  of  the  reservoirs.   The 
first  condition  allows  the  steady  state  to  be  approached.   The  sec- 
ond condition  prevents  energy  from  being  reflected  at  the  ends  of 
the  reservoirs,  travelling  in  the  opposite  direction,  and  causing 
the  entire  system  to  approach  thermal  equilibrium.   As  a  rough  esti- 
mate, we  may  in  general  choose  T  such  that 


^0  «  T  «  ^z  II--^ 


where  v  is  the  phonon  speed  and  L^  is  the  total  length  of  the  crys- 
tal along  the  z  axis.   Since  L^  is  to  be  made  infinite,  T  may  be 

chosen  as  infinite. 

Once  p(t)  has  been  found  the  statistical  average  of  the 
steady  state  value  of  the  current  <J >  may  be  found  by  evaluating 


<J>=  Tr  J  p(t)  II-5 

where  J  is  the  quantum  mechanical  current  operator.   The  actual  deri- 
vation of  p(t)  and<J>is  lengthy  and  quite  difficult.   Since  this 
derivation  is  adequately  discussed  in  references  noted  previously, 
only  the  result  will  be  given  here. 

For  this  model,  the  average  value  of  the  heat  current  is 

,26 

given  in  the  Kubo  formalism  by 
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<X>=  —  /  ds   /  dX  Tr  [J(0)  J   (t+iXTi)p-]  AT  .   II-6 
0     0  ^         U    z 

In  this  equation,  A  represents  the  cross-sectional  area  of  the  sample, 
^■^z  ^^^   temperature  difference  between  the  reservoirs,  T  the  average 
temperature,  and  X  is-  a  parameter  having  units  (energy)"  .   p   is  the 
equilibrium  density  matrix  given  explicitly  by 

-> 

and  J  is  the  local  heat  flux  operator  averaged  over  the  sample  of 

thickness  L_  (to  be  discussed  in  the  next  section).   In  the  Heisenberg 
picture,  the  time  dependence  of  the  operator  J  (t+iXTi)  is  given  by 


J  (t+iX-h)  =  exp(iHt/h)  exp(-XH)  J  (0)  exp(XH)  exp  (-iHt/h).   II-8 


If  one  substitutes  Eqs.  II-7  and  II-8  in  Eq.  II-6  and 
evaluates  the  trace  using  Dirac  bra  and  ket  notation,  one  finds 


<J>  =  ^  /  ds   /  dX  I   <5,|j|n>  <nlj  |£>  exp[-X(E  -Ej]  exp[is(E  -E,)/n] 
0     0    £,n  "   ^  "  ^ 


exp(-B  E  ) 

^L_  AT  II-9 

I   exp(-6  E,  )      ^ 


where  J  =  J(0).   |£>  and  |n>  represent  eigenvectors  of  the  Hamiltonian 

having  energy  eigenvalues  E  and  E,  and  satisfy  the  orthonormality 

11       x* 
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condition 


<Mn>  =  6^^        .  II-IO 


The   integrals   over   s   and   X   can  easily  be  performed    to   yield 


<J>  =     -~       I     <£ljln>  <n|j    \Z>   6/\   h\   p..   AT        11-11 


where 


exp(-6E^) 
P0£  =      I   exp(-6K  )        •  11-^2 

k  '^ 


Eq.  11-11  is  the  specific  equation  which  will  be  used  in  the  evalua- 
tion of  the  heat  current  and  consequently  the  thermal  conductivity. 

2.   The  Average  Heat  Flux  Operator 

28 
Hardy   has  shown  that  the  local  heat  flux  operator  can 

be  obtained  by  requiring  that  it  satisfy  the  continuity  equation 


H(x)  +  V'J(x)  =  0   .  11-13 


H(x)  is  the' time  derivative  of  the  Harailtonian  density  operator 
which  is  related  to  the  Hamiltonian  H  of  the  system  through  the 
Heisenberg   equation  of  motion 


H(x)    =    (ih)"^    [H(x),H]    .  11-14 
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For  a  general  type  Hamiltonian  of  the  form 

2 


i      1 


Hardy  has  shown  that  the  harmonic  contribution  to  the  local  heat 
flux  operator  is  given  by 

^(x  ) 
J  (x)  =  I  I      A(x-x^)(x^-x  )(iti)~l  [    ^   ,V(x  )]   11-16 
k,j  -■         He       -' 

where  the  sum  over  k  and  j  is  to  be  performed  over  the  entire  crys- 

-> 
tal.    In  these  expressions,  x,  is  the  position  vector  to  the  equi- 

th  "*■  "*" 

librium  position  or  lattice  site  of  the  k   particle,  p(x  )  its 

->• 
momentum  operator,  and  m,  its  mass.   V(x.)  is  the  potential  associ- 

k  J         '^ 

ated  with  the  j    atom  and  in  the  harmonic  approximation  is  given 
by 


^(V  =  i  K  I 'as  ^Vv  \^v  "3^v      "-^^ 

-"      a,3  Ji      -"        -■ 

where  $  „(x.)  is  the  so-called  force  constant  tensor  obtained  by 
a3  J 

expanding  the  potential  in  a  Taylor's  series  about  its  equilibrium 
position  and  retaining  only  the  quadratic  term.   Quite  generally, 

-*■  -y 

^  o(x.)  is  a  symmetric  function  of  x.and  is  also  symmetric  with 
a3  J  J 

->-  -> 
respect  to  interchange  of  a  and  3.   The  vector  U(x.)  represents  the 

displacement  of  particle  j  from  its  equilibrium  position.   Greek 

subscripts  are  used  to  denote  Cartesian  components  of  vector  and 

tensor  quantities  throughout.   Finally,  A(x-x,  )  is  a  coarse-grained 
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delta  function  used  by  Hardy  which  he  has  chosen  to  specifically 
represent  by  the  formula 


,->->■  1 2   2 


A(x-Xj^)  =  ^3/2  p   exp[-|x-x^|  /£^]  .        11-18 


In  Eq.  11-18,  Z   is  distance  large  microscopically  but  small  macro- 

scopically.   The  function  A(x-x, )  is  normalized  such  that 

k 


/dx  A(x-Xj^)  =  1   .  11-19 


By  using  the  commutation  relation  for  the  position  and 
momentum  operators  and  Eq.  11-17,  one  can  show  that  for  i  j"^  j 


[p^(x,),V(x.)]  =-ih  I   '^.^gCx.-x.)  pjx.)  UpCx.)  .     11-20 


Hence  Eq.  11-16  becomes  ■ 


.   11-21 


For  further  development  of  this  equation,  it  is  necessary  to  trans- 
form  the  operators  p(Xj^)  and  U(x.)  to  the  normal  mode  or  phonon 
representation.   For  the  case  of  a  single  atom  per  unit  cell,  such 

a  transformation  is  accomplished  by  the  well-known  transformation 

,.29 
equations 


U(x,)  =   l^  ^,     [  B^^xJ  a(q,s)  +  t^^a:>    ,Uts)] 


J,s  V  2m^  (^(q.s) 
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and 


11-22 


->■  -> 


i,hoo(q,s) 
P(x^)    =  -i     J\/-^ [   S^^xp    a(J,s)    -  r^«(J^)    a+(J,s)]     . 

q,s 


Here  aj(q,s)    is    the   frequency  of  phonons   having  wave  vector   q   and 

t  ->■ 
polarization  s.      a    (q,s)    is    the   creation  operator   for    these  phonons 

and   a(q,s)    the   corresponding  destruction  operator.      These   operators 

satisfy    the   commutation  relation 


[a(q,s),a    (q',s')]    =   6-^,    6       ,         .  11-23 

qq         ss 


Virtien  operating   on   the   state  having   N   phonons  of  wave  vector  q   and 
polarization  s,    they  have   the   effect 


I    (q,s)  |n(J,s)>   =    ^(J,s)4-1    I    N(q,s)+1> 


and  11-24 


.(J,s)  |N(q,s)>   =    A(q,s)     1    N(J,s)-l> 


where  |N(q,s)+l>  symbolically  means  the  state  in  which  N(q,s)  has 

t  -^ 
been  increased  by  unity  and  so  on.   Additionally,  a  (q,s)  is  the 

Herraitian  conjugate  of  a(q,s)  and  B  ^  (-<q)    is  the  complex  conjugate 

of  B   (x„).   Finally  the  B   (x„)  satisfy  the  equation 

3,k   i^ra^n^ 
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and  are  orthonormal  in  the  sense  that 


->* 


B  ^^(x,)  .  B^  ^  (xj  =   6^,    6      ,   .  11-26 

^  X,     qq'   ss' 


If   Eqs.    11-22  are  used   in   Eq.    11-21,   one  obtains 

->     ->- 

'     k.j      a  3     s,s'    j^j,  ^       "^     ^    ■    y^ 

11-27 

?'^'     ^  4s    ^  ^  X    ^  Vs'     ^  Js    _^  .     _^ 

"^    ^\         ^V^3       (^j)^(q'.s')a'(q,s)-B^  (Xj^)Bg    (xja'(q',s')a(q,s)] 


In  obtaining  Eq.  11-27,  those  terms  which  do  not  conserve  energy 
have  been  neglected.   As  evidenced  by  Eq.  11-11,  they  make  no  con- 
tribution  to  <J>. 

The  average  value  of  the  heat  flux  operator  is  given  by 

■*■     1   f  ->  ->■  ,->■ 

■^2  "  v"  J  '^^  ^^2^^^  11-28 

s 

where  the  integral  is  to  be  done  only  over  the  sample  and  V  repre- 
sents the  volume  of  the  sample.   For  those  atoms  lying  outside  the 
sample,  the  coarse-grained  delta  function  causes  the  integrand  to 
become  negligibly  small.   Furthermore,  if  one  assumes  that 
*^a3  k~^i^  is  negligibly  small  when  x  -x .  is  larger  than  a  few 
lattice  spaces,  the  average  heat  flux  operator  can  be  written  with 
negligible  error  as 


^2  =  iv~    r    ^  .  ^ss'^'^'^  a'^cj.s)  a(J',s')     11-29 


s   -^  ^,    s,s' 

q,q 
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where 


^ss'(^''^'>=  I.   I  (-r-j)   f-l  ^  \        (-,)  B3   (X.).   11-30 

The  notation  I      is  used  to  indicate  that  the  sum  is  to  be  done  only 
over  the  sample. 

Before  concluding  this  section,  it  is  necessary  to  note 
some  properties  of  the  quantity  G   j(q,q')  which  will  be  used  in 
future  calculations.   First,  inspection  of  Eq.  11-30  reveals  that 


Ggg,(q,q')  =  -G*,g(q',q)   .  11-31 


Second,  if  one  performs  the  sum  over  j  in  Eq .  11-30,  Eq.  11-25  may 
be  used  to  prove  that  the  z  component  of  the  resulting  expression 
is  independent  of  z„   provided  a)(q,s)  =  a)(q',s),  which  is  the  only 
case  of  interest.   This  rather  intuitively  obvious  result  just  means 
that  there  exist  no  sources  or  sinks  for  the  current  along  the  z 
axis. 

3.   Equations  of  Motion 

The  final  bit  of  formalism  to  be  discussed  is  the  difficult 

-> 

q  g  "> 

problem  of  determining  the  quantities  B   (x„)  which  appear  in 

30 
Eq.  11-30.   The  method  to  be  used  is  due  to  Maradudin   and  employs 

the  Green's  function  technique. 

By  defining 

->         1/2  -> 
B^^(x.)  =  m.    C^^(x.)  A(J,s)  11-32 
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and  substituting  in  Eq.  11-25,  one  finds 


qs  qs 


2,-* 


which  is  just  the  equation  of  motion  of  the  i       atom  of  the  crystal. 
Physically,  C^^Cx.)  is  the  displacement  of  the  I        atom  from  its  equi- 
librium  position  provided  only  the  q   mode  is  excited. 

If  we  recall  that  the  crystal  under  investigation  has  host 
atoms  of  mass  M  as  well  as  impurities  of  mass  M'  within  the  sample, 
we  may  re\rc±te   Eq.  11-33  in  the  form 


k,p  k,p 


or   in   the  more   abbreviated  matrix  notation  as 


L  C"^^  =   6   L  C'^^  11-35 


where 


^aB^^rV   =  ^  "^  ^a^  h\.  -  '^aB^VV  ""^^ 


and 


^   W^iA^   =   ^  ^  "'   I  ^aB  \n     \n        "  ""'^ 


In   these   equations,    e   is    the  mass  defect  parameter   given  by 

e  =   1  -   M'/M  and   Xj^      is   used    to   denote    the  position   of    an   impurity, 


23 


It  may  be  noted  that  6L  and  thus  the  right-hand  side  of  Eq.  11-35 
are  zero  if  no  impurities  are  present,  and  one  just  obtains  the 
equation  of  motion  of  the  isotopically  pure  crystal.   For  the  case 
in  which  impurities  are  present,  6L  is  non-zero  only  when  k  and  i 
denote  positions  of  impurities. 

Maradudin  has  shown  that  the  solution  to  Eq.  11-34  can  be 
written  in  the  form  of  a  plane  wave  incident  on  the  sample  plus  a 
wave  scattered  by  the  sample,  i.e. 


->         qs        ^ 
C^^(x^)  =  Cq  (t^)   +   W^^(Jj^)  11-38 


->  *  -^Xi  s 

where  C   is  the  incident  plane  wave  and  W   the  scattered  wave. 
The  incident  wave  C_   is  a  solution  to  Eq .  11-34  when  no  impurities 


are  present  and  thus 


qs 

L  Cq  =  0   .  11-39 


If    Eq.    11-38    is    substituted    in   Eq.    11-35,    one    finds 


->■  •       qs  ^ 

L  w"^^  =  5  L  Cq     +  5  L  w'^^      .  11-40 


Finally  if  the  matrix  G  =  L   is  defined,  Eq.  11-40  can  be  written 
in  the  form 


->        qs       -> 
W''^  =  G6LCq  +  G6Lw'^^   .  11-41 
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The  matrix  G   is   known  as    the  Green's   function  and   has  been   shown  by 

30 
Maradudin   to   be  given  by 


^a3^''£'V  ~ 


M(2tt) 


3   /dq 


q^q 


a^3 


2      2   2 

.  (jj  -V  q   +i6 

Li 


2      2    2^.. 


+  6 


aB        2     2   2^.^ 
CO  -v^q  +16 


exp [iq*x(£,k) ] 


11-42 


where   "a"   is    the   lattice  constant,   v      and  v      the  group   velocities   of 
phonons   in   the   longitudinal   and    transverse  modes,    6   a   small  real 
number,    and 


x(5,,k)    =  x^  -  Xj^ 


11-43 


In  writing  Eq.  11-42,  it  has  been  assumed  that  the  acoustic  approxi- 
mation is  valid,  i.e. 


a)(q,s)  =  V  q   . 
s 


11-44 


The  integral  in  Eq.  11-42  can  be  performed  by  integration  in  the 
complex  plane.   The  resulting  expression,  however,  is  quite  complicat- 
ed and  will  not  be  given  here.   Since  our  ultimate  concern  will  be 
for  the  case  in  which  v  =  v  ,  we  note  that  for  this  case  Eq.  11-42 

J.  Li 

becomes    the   rather   simple   expression 


'"s'^^'^j'  '  ;^ 


_     exp[icox(£,  j)/v]    . 
2  x(£,j)  aB 


11-45 
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where  v  =  v  =  v  .   It  should  be  emphasized,  however,  that  these 
velocities  have  been  set  equal  merely  to  simplify  calculations  and 
that  there  is  no  fundamental  difficulty  in  using  the  more  general 

Green's  function.   In  fact,  calculations  have  been  performed  by 

31 
Allen  and  Powell   for  the  case  where  v  ^   v  ,  though  for  a  simpler 

Li         X 

problem, 

-> 

"^  s 

The  ultimate  aim  is  to  solve  Eq.  11-40  for  W^  .   A  formal 
solution  can  be  obtained  by  writing 


->  _      qs 

W'l^  =  [1-G6L]    G5LCq   .  11-46 


For  purposes  of  calculation  we  will  expand  the  matrix  [1-G5L] 
which  yields 

->  ->• 

->        qs  qs 

W^^  =  G6LC   +  G5LG6LC   +  ...  11-47 


The  first  term  in  Eq.  11-47  represents  the  first  Born  approximation 

and  subsequent  terms  are  of  higher  order  in  the  Born  series. 

-> 
"*XIS  "*■ 
Once  the  value  of  the  scattered  wave  W^  (Xp)  has  been 

-> 
found,  the  quantities  B^^(xp)  can  be  immediately  obtained  as  can  be 

seen  from  Eqs.  11-32  and  11-38.   Quite  formally,  one  may  write 


B^^x.)  =m.   A(J,s)  [C   (x.)  +W^^x  )]         11-48 

-y 
where  W^^(x.)  is  given  by  Eq.  H-46  or  11-47.   The  normalization 

constant  A(q,s)  can  be  determined  by  use  of  Eq.  11-26. 
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This  section  completes  the  formalism  required  for  a 

-)- 
complete  calculation  of  <J>  given  by  Eq.  11-11. 

4.   Synthesis  of  Results 

The  results  of  the  preceding  three  sections  may  be  sum- 
marized by  noting  that  we  have  presented  formal  methods  for  calcu- 
lating all  quantities  necessary  to  evaluate  the  steady  state  aver- 
age  heat  current  <J>.   Specifically,  in  section  two  the  average  heat 

flux  operator  was  obtained  and  in  section  three  a  method  was  given 

->■ 
->-qg  ->- 
for  finding  the  quantities  B^  (x,)  upon  which  that  operator  depends. 

In  this  section,  the  results  will  be  combined  as  much  as  possible, 

without  loss  of  generality,  in  order  that  an  expression  for  <J>  which 

is  perhaps  more  amenable  to  calculations  may  be  obtained. 

Recalling  that  our  interest  is  only  in  the  z  component  of 

<J>,  we  have  upon  substitution  of  Eq .  11-29  into  Eq.  11-11 


<-J,>  =  _i[AhL  I  I  I     G^,3.  (?.?')  G(q' ',?"•) 

HKL    V   q,q         , ,   , , , 

q",q'" 


6(0j(q,s)-Oj(q',s'))<£la(q',s')a^(q,s)|n> 


X   <n|a(q'",s'")a'^(J",s")|£>  p^„AT   . 


.  ^^^J^t   „Nl„^  11-49 


or 


z 


If  we  define 


"ss's"s'"^^'^''^"'^'"^  =   I  <S.|a(q',s')a'^(q,s)|n> 

£,n 


X  <n|a(q'",s"')a^(q",s")|£>  p^^   , 


11-50 
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Eqs.  11-10  and  11-24  may  be  used  to  show  that 


M^,i,ri^,,,(q,q',q",q'")  =  <5,„i  1 1^^,^,  ,6^, ,  ,5;^,;>, , 


e 


ss's"s'"'^'^  '^   '^    '         "ss"'  s's"  qq'"  q'q"  ,  ^  ,,2    ' 

(e  -1; 


->.  ->-. 


q  i    q',  q"  i    q"',  s  j^  s',  s"  4    s'" 


M3^,3,,3,,,(q,J',J",J'")  =    ,   ^   g, +  — I 2  ^?^"^ss" 

^   ^  (e^-l)(e^  -1)    (e^-1)^  ^^        ^^ 

q  =  q',  q"  =  q'",  s  =  s',  s"  =  s'" 


Mgg,g,,g,,,(q,q',q",q'")    =   0      ,    otherwise  11-51 


where 


^    -  kT  •  ^^    ^^ 


In  deriving  Eqs,  11-51,  use  has  been  made  of  the  fact  that  since 
phonons  obey  Bose  statistics 


I   N^(q,s)  Pqj^  =  <N(J,s)>  =     1 
^  e^-  1 

and  11-53 


2  ^ 

^  N  (J,s)  p    =  <N^q,s)>  =  ^\^      . 


Use  of    Eqs.    11-50  and    11-51   in   Eq.    11-49   yields 
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<J  > 
z 


11  Ah 

q,q 


11-54 


6(to(q,s)  -  aj(q',s'))  AT 


where  use  has  been  made  of  Eq.  11-31. 

If  one  wishes,  the  explicit  expression  for  G   i(q,q')  given 


ss 


^s 


by  Eq.  11-30  and  for  B^  (x.)  given  by  Eqs.  11-48  and  11-47  may  be  sub- 
stituted into  Eq.  11-54,   No  further  simplification  results,  however, 
in  the  most  general  case  and  this  will  not  be  done  here.   Eq.  11-54 
will  be  used  in  Chapter  IV  where  explicit  calculations  of  the  current 
are  undertaken. 


CHAPTER   III 
THE  SCATTERiriG  PROBLEM 

The  purpose  of  this  chapter  will  be  to  obtain  a  solution 

-)- 

for  the  scattered  wave  W   (Xj,)  discussed  in  the  previous  chapter  and 
the  B^  (xp)  appearing  in  Eq.  11-30.   It  will  be  assumed  that  the 
impurities  are  uniformly  distributed  throughout  the  sample.   It  is 
further  assumed  that  the  wavelengths  of  the  phonons  are  much  longer 
than  the  distance  between  successive  impurities.   The  significance 
of  the  long  wavelength  assumption  in  the  calculation  of  the  current 

will  be  discussed  in  the  following  chapter.   In  addition,  the  solu- 

-y 

tions  for  the  b'^  (xp)  are  found  only  outside  the  scattering  region, 
but  it  will  later  be  ai'gued  that  these  solutions  are  sufficient  for 
calculating  the  current. 

1.   Scattering  by  a  Single  Plane  of  In-purities 

We  begin  by  considering  a  single  plane  of  atoms  located 
in  the  infinite  plane  z  =  0.  A  certain  fraction  O  of  these  atoms 
ai'e  impurities  and  these  are  uniformly  distributed  among  the  host 
atoms .  A  plane  wave  of  the  form 


->^s  -^  -^  -^       -y  ^ 

Cq  (x^)  =  exp(iq'x^)  A(q,s)  III-l 


where  q  is  the  wave  vector  and  X(q,s)  the  polarization  vector,  is 
incident  from  the  left  on  the  plane.   Since  the  impurity  plane  will 
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be  treated  in  the  continuum  limit,  there  is  perfect  azimuthal  sym- 
metry about  an  axis  normal  to  that  plane.   Thus  the  propagation  vector 
of  the  incident  wave  can  be  chosen  to  be  in  the  yz  plane  and  make  an 
angle  a  with  the  z  axis.   (See  fig.  2.) 

If  we  denote  the  first  term  in  the  expansion  of  Eq.  II-U6 
by  W-,  jthen 


^  ->      .    ^,-> 


qs  2  3 

^1  ^^^^  =  r-2 

r 


^'^^  f^    ^    -ceo  a^   r   exp[iqx(£n  )]  exp(iq'Xn^)  X(q,s)    III-2 
W.  (xj  =  —^ —  I     r r . 


^^     ^  xiln    ) 


In  Eq.  III-2  explicit  use  has  been  made  of  Eqs .  II-3T,  11-!+^  and 
II-i45.   If  it  is  further  assumed  that  the  distance  between  succes- 
sive impurities  is  much  smaller  than  the  phonon  wavelength,  the 
sum  over  r  in  Eq.  III-2  can  be  replaced  by  an  integral  over  the  im- 
purity plane  with  the  result 


qs  2  exp[iq(x(5,n  )+yn  sin  ct)] 

t    it)   =     --^^  I   dx^   -^-^ '—^ X(q^s).  III-3 

^   ^      Utiv^      ^       x(£n  ) 


In  Appendix  A,  it  is  shown  that  integrals  of  this  form  can  be  approxi- 
mated by  the  following  formula  which  is  a  consequence  of  the  method 
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of  stationary  phase   : 


//g(x',y')  exp[iqf(x',y')]dx'  dy'  ^ 

^  III-U 

a . 

~  I     I  W=-     sU.^y.)    exp[iqf(x  y  )1  . 
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Figure  2 

Orientation  of  vectors.   The  incident  plane  wave  denoted 
by  C_  is  chosen  to  lie  in  the  yz  plane  and  make  an  angle  a  with  the 
z  axis.   The  impurity  plane  is  the  z  =  0  plane,  x^   and  x^^  repre- 
sent vectors  from  the  origin  to  the  Ip.-.tice  sites  of  a  host  atom 
and  an  impurity,  respectively. 
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In  Eq.  III-U,  X.  ■and   y.  are  points  in  the  domain  D  at  which  f(x',y'; 
is  stationary  and 


a.  = 


9!l 

9x'^ 


t   I    <J 


3!l 


3"f 

,   ,  '  ^j    3x'3y' 

X.  ,y . 
J  J 


,  ,  .  III-5 

X.  ,y . 
J  J 


The  number  O.    is  given  by 


a.  =  1 
J 


a.  =  -1 
J 


,  for  a. 6.  > 
J  J 


Y 


a.  >  0 


,  for  a.6.>  Y.  ,a.  <0 
J  J     J    J 


III-6 


a.  =  -1 


,  for  a.3.  <  Y.   . 


If  Eq.  Ill-';  is  appliei  to  Eq.  III-3,  one  finds  after  a  tedious  but 
straightforward  calculation 


as 


-  ietjqa 


->-+-»■   V   ->, 


\   ^H^  =   2lcos  a|  ^^P(i^i-^£)  ^^1'^)  • 


III-7 


The  upper  and  lover  signs  on  the  vector  q  are  to  be  taken  when 

z.  cos  a  <  0  and  z„  cos  a  >  0,  respectively,  and  qT   and  q  are  given 


i3y 


q     =  q(sin  aj    -  cos  otk) 
q  =  q(sin  aj    +   cos   ok) 


III-8 


Throughout  the  remainder  of  the  discussions  q,  will  be  used  to  denote 
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the  vector  q  i  +  q  j  -  q  k  where  q  ,  q^  ,  and  q  are  the  components 
of  the  vector  q. 

A  thorough  physical  discussion  of  the  scattered  wave  will 
be  given  later.   For  the  moment,  however,  it  is  sufficient  to  note 
that  a  plane  wave  incident  on  the  impurity  plane  will  produce  a  re- 
flected plane  wave  at  a  reflection  angle  equal  to  the  incidence  angle 
and  a  plane  wave  scattered  in  the  forward  direction  without  deviation. 
This  can  be  seen  by  examination  of  Eq.  III-T.   Clearly  the  case 
z.  cos  a  <  0  represents  the  reflected  wave  and  z.  cos  a  >  0  the  for- 
ward scattered  wave. 

The  first  Born  amplitude,  represented  by  the  coefficient 
of  the  exponential  in  Eq.  III-T  should  be  a  good  representation  of 


the  actual  amplitude  for  small  values  of  q,  e,  a   and  for  a  not  too 
2' 


TT 

close  to  — .   It  may  be  noted,  however,  that  this  amplitude  is  singu- 


lar  at  a  =  —  .   The  soxirce  of  this  singularity  can  be  traced  to  the 
fact  that  as  a  ->  —  ,  the  scattering  becomes  very  strong  and  the  Born 
approximation  is  no  longer  valid.   Since  ultimately  the  intensity 
of  this  scattered  wave  will  have  to  be  integrated  over  all  values  q, 
the  singularity  must  be  removed.   The  removal  of  this  singularity  can 
be  most  effectively  accomplished  by  consideration  of  the  higher  order 
terms  in  the  Born  expansion.   For  example,  the  second  Born  term  in 
Eq.  II-U7  can  be  written 


x(£n  ; 
r 


3U 


and  if  one  assumes  that  the  first  Born  term  given  by  Eq.  III-T  repre- 
sents an  adequate  solution  to  the  amplitude  inside  the  scattering 
region,  a  calculation  identical  to  the  calculation  of  Eq.  III-7  can 
be  performed  with  the  result 


2  cos  a 


exp(iq^«Xj^)  A(q,s)  .         III-IO 


III-ll 


Clearly  this  process  can  be  continued  to  produce 


-*■  qs  7 

W^  (x^)  =  I   W^  (x^)  =  expdq^'x^)  X(q,s) 


X   [-ix(q)  +  (-.ix(q))^  +(-ix(q))^  +  .  .  .] 


where 


x(q)  =  o-r^^S-   •  "1-12 

^    2|cos  a| 


If  we  now  assume  that  this  expression  can  be  resummed  to  produce  a 
resiilt  in  the  form  of  Eq.  II-U6,  we  have 

and  from  Eqs ,  III-l  and  11-38, 

C^^(J^)  =  exp(i^^^)  ta,s)   -   -i^  exp(i^.x^)  A(J,s)  .  III-H; 

l+ix(q) 
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Eq.  III-13  will  be  used  to  represent  the  value  of  the  wave 
scattered  "by  the  single  plane  of  impurities.   Admittedly,  one  has  no 
rigorous  assurance  that  Eq.  III-7  represents  a  valid  solution  for  the 
scattered  wave  at  points  lying  in  the  impurity  plane  as  is  discussed 
in  Appendix  A  and  thus  the  above  method  of  calculating  the  higher 
order  terms  in  the  Born  series  is  perhaps  questionable.   However, 
Eqs.  III-13  and  III-lU  do  possess  the  following  meritorious  proper- 
ties: 

(i)    They  contain  no  singularities. 
(ii)   They  reduce  to  the  first  Born  approximation 
for  small  values  of  x(q). 

TT 

(iii)   They  predict  that  as  a  ->■  j  ,  the  intensity 
of  the  transmitted  wave  vanishes  and  the 
intensity  of  the  reflected  wave  approaches 
unity.   This  is  analogous  to  the  well-known 
optical  case  of  electromagnetic  waves  re- 
flected and  transmitted  by  a  dielectric 
interface. 

(iv)   Finally, they  predict  values  of  the  B  (x.) 
which  satisfy  the  orthogonality  condition 
given  by  Eq.  11-26. 

In  short,  though  the  technique  used  to  calculate  these  higher  order 
terms  lacks  mathematical  rigor,  it  does  produce  results  which  appar- 
ently satisfy  all  mathematical  and  physical  requirements.   In  addi- 
tion, for  values  of  the  temperature  which  will  be  considered,  x(q) 
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will  be  very  small  except  for  values  of  a  exceedingly  close  to  —  . 

For  these  values  it  will  be  seen  that  the  current  is  approaching  zero. 

-> 

We  therefore  conclude  that  use  of  Eq.  III-14  for  the  C   (x.)  cannot 
cause  significant  error. 


2.   Scattering  by  Many  Planes 

In  Eq.  1X1-1^+,  we  have  essentially  found  the  reflection 
and  transmission  coefficients  for  a  single  plane  of  impurities.   De- 
noting these  by  t  and  r,  we  have 


t  = 


r  = 


l+ix(q) 

III-15 

l+ix(q) 


~^s  / ~^  \ 
In  this  section  we  find  the  values  of  C "  (x. )  and  consequently 

B   (x.)  for  the  case  in  which  scattering  can  occur  from  N  planes. 
The  planes  are  uniformly  spaced  a  distance  "a"  apart  along  the  z 
axis  in  the  length  L  of  the  sample.   Multiple  scattering  is  account- 
ed foro   Though  several  modifications  have  been  introduced,  our  ana- 

33       3^  35 

lysis  will  be  similar  to  that  of  James,   Prins,    and  Darwin   who 

obtained  formal  equations  for  the  total  electromagnetic  wave  reflect- 
ed from  and  transmitted  by  any  given  plane  of  atoms  in  a  lattice  for 

the  case  of  Bragg  reflection.   The  solution  of  these  enuations  is 

?6 
similar  to  that  obtained  by  Mooney   for  the  case  of  optical  filters. 

We  begin  by  numbering  the  planes  consecutively  along  the 

z  axis  beginning  with  the  first  plane  at  z  =  0  and  ending  with  the 
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N   plane  at  z  =  (N-l)a.   (See  fig.  3.)   V7e  assume  that  a  plane  wave 
having  unit  amplitude  and  unit  phase  at  z  =  0  is  incident  from  the 
left  (i.e.  cos  a  >  O).   The  amplitude  and  phase  of  the  total  wave  re- 
fleeted  from  and  of  the  wave  transmitted  down  to  the  k   plane  are 
denoted  by  R  and  T  ,  respectively.   These  values  refer  to  the  ampli- 
tude and  phase  just  to  the  left  of  the  k   plane.   Except  for  the 
first  plane  and  the  N   plane,  T  consists  of  that  pai't  of  T 
transmitted  through  the  k-1   plane  as  well  as  that  part  of  R  which 

K, 

is  reflected  from  the  left-hand  side  of  the  k-1   plane.   Similarly, 
R,  consists  of  that  part  of  T  reflected  from  the  k   plane  as  veil 
as  that  part  of  R-,  .-,  which  is  transmitted  through  the  k   plane. 
The  contributions  to  the  waves  for  k  =  1  and  k  =  N  can  be  determined 
by  similar  arguments.   Thus  we  have 


T^  =  l 


R  =  T  r  +  R  te  ^ 


III-16 


T,  =  T,   te^*  +  R,  re^^* 
k    k-1       k 


for  k  =  2,3...,N-1 


icf) 
k'  ■  ^'k+1'" 


R^  =  T,_r  +  R,_.,te 


ict       2x6 
^N^^N-l^^^V^ 


R  =  rT 
N     N 
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N 


Rk 


z=o 


Tk 


Z=(k-i)a 


Z=([M-I)a 


Figure  3 

Scattering  by  many  planes.   Planes  are  n-umbered  from  left 
to  right  beginning  at  the  origin,   C_  represents  a  plane  wave  inci- 
dent from  the  left  makirig  an  angle  a  with  the  z  axis.   R,  and  T, 
denote  amplitudes  and  phases  of  the  total  wave  reflected  from  and 
transmitted- down  to  the  k   plane.   That  is,  they  refer  to  values 
attained  just  to  the  left  of  the  k^  plane.   Only  a  cross-sectional 
view  of  the  planes  has  been  depicted. 
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where  <p  =   aq|cos  a  |  is  a  phase  factor  which  accounts  for  the  dif- 
ference in  path  travelled  by  the  various  contributions  to  the  waves, 
E^s .  III-I6  can  be  solved  for  R  and  T  in  terms  of  R 


and  T  ^  and  written  in  the  form 
k-1 


T^  =  l 


R    -  T    r 
„  _   k-1    k-1 

^'  "IT* 


\         r  t     \-l         t   V] 


for  k  =  2,3...N 


R   =  rT 

N     N 


or  in  matrix  notation  as 


T^  =  l 


III-IT 


'\.] 


k+1' 


R, 


=  M 


ho 


for   k  =  1,2...N-1 


R  =  rT 
N     N 


In  Eq.  III-I8,  M  is  the  two-by-two  matrix 


III-I8 


(l+ix)e-^* 


M  = 


-ixe 


i4> 


ixe 


-i4) 


(l-ix)e 


icf. 


III-I9 


In  obtaining  Eq.  III-I9,  Eq.  III-I5  has  been  used.   Finally  Eq. 
III-I8  can  be  applied  successively  to  obtain  the  following  equat 


lons! 


T^=l 


,B..S 


R, 


^^=M^-^"^ 


^n/ 


III-20 


R  =  rT 
N     N 


Ui 


N-1 
We  denote  the  elements  of  the  matrix  M    by 


M^-l  = 


a  h' 
c   d 


III-21 


and  have  from  Eq.  I I I- 20, 


T^  =  l 


R  =  rT 


^N  =  ^1  ^  bT^ 


Tjj  =  cR^  +  dT^      .  III-22 


These  equations  can  be  easily  solved  for  R  and  TL,   yielding 


R,  =  ^^1^^ 
1    a  -  cr 


N     a  -  cr 


R  and  T  represent  the  values  of  the  amplitude  and  phase  of  the 
total  wave  (including  all  posible  multiple  reflections)  reflected 
from  the  first  plane  and  transmitted  down  to  the  N   plane. 

For  further  Sxmplification,  we  note  that  the  matrix  M  of 
Eq,  III-18  has  a  unit  determinant.  For  such  matrices  it  has  been 
shown   that  the  elements  of  the  matrix  M  can  be  written  in  terms 


h2 


of  the  elements  of  the  matrix  M  and  the  Chebyshev  polynomials  of  the 
second  kind.   Specifically, 


M^  = 


1   III-24 
i^^21  Vl^Q)  ■  ^22Vi^^)-"n-2(^) 


where  the  M.   are  the  elements  of  the  matrix  M  and  0  is  given  by 


0  =  1/2  (M^^  +  M^g)    .  III-25 


The  U^_^(0)  are  Chebyshev  polynomials  of  the  second  kind  and  are 
given  explicitly  by  the  formula 

U^_^(0)  .   sin  [(N.1)  cos-^  01  ^  ^^^_^^ 


VTTe^ 


These  polynomials  fui'ther  satisfy  the  recursion  relation 


U  (0)  =  2  U    (0)  -  U.  .(0)  ■  III-2T 


and  the  orthonormality  condition 


/  u  f0)  u  (0)  Vi  -  0^  -      ^ 


_^  'M^-'  ^N^-'  ^  -  -  -  ^®  =  2  hm    '  ^"-28 


Applying  these  results  to  find  the  elements  of  the  matrix 


1^3 


given  by  III-21,  we  have 


a  =  (l+ix)  e"^*  U^_2(9)  "  "n-3^®^ 
b  =  ix  e"^*  U^_2(e) 


c  =  -ix  e^*  Uj^_2(G) 


d  =  (1-ix)  e^*  V2^®^  ~  Vs^®^        ■  ^"-^^ 


where 


0  =  cos  (j)  +  X  sin  (j)    .  III-30 


Finally,  if  we  denote  the  values  of  the  amplitude  and  phase  of  the 
total  wave  transmitted  by  the  N  planes  by  T(q)  =  ^T  ,  we  have 


T(q)   = :—      .  III-31 

(1+ix)  U^j_^(0)  -  e^*  U^_2(0) 


Similarly  the  total  wave  R(q)  =  R  reflected  by  the  N  planes  is  given 
by        '      • 


R(q)   =  ^-^ — .■  III-32 

(l+ix)U^_^(0)  -  e^*  U^_2(0) 
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In  obtaining  Eqs .  III-31  and  III-32,  Eqs .  III-23,  III-2T  and  III-29 
have  been  used. 

The  analysis  above  has  been  for  the  case  in  which  the  in- 
coming wave  was  incident  from  the  left  or  for  cos  a  >  0.   However, 
since  T(q)  and  R(q)  depend  only  on  the  absolute  value  of  cos  a  , 
Eqs.  III-31  and  III-32  are  also  valid  for  cos  a  <  0.   We  therefore 
conclude  that  a  plane  wave  having  unit  amplitude  and  phase  incident 
from  the  right  onto  the  N   plane  produces  a  wave  transmitted 
through  thf;  first  plane  whose  amplitude  and  phase  are  given  by 
Eq.  III-31,  and  a  wave  reflected  from  the  N   plane  whose  amplitude 
and  phase  are  given  by  Eq.  III-32. 

We  now  wish  to  apply  the  above  results  to  find  the  ampli- 
tude  and  phase  at  any  arbitrary  point  x„.   Using  the  notation  of  the 
previous  section  and  recalling  that  part  of  the  phase  of  the  total 
wave  is  accounted  for  by  the  factor  exp(iq  •  x  ),  we  have 

->■ 
C^  (Xp)  =  T(q)exp[-iq(N-l)alcos  a|]exp(iq*x  )  A(q,s) 


for   z.cos  a  >  0 


;^s 


-)-->.->,->-  ,     ,->-^     ,->•->■»->,-> 


C^  (x^)  =  exp(iq'X^)  X(q,s)  +  R(q)  exp(iq»x^)  X{q,s) 


;^s 


for  Zj,  <  0  and  cos  a  >  0 


->-  ->■  N  :>,-> 


C^  (Xj^)  =  exp(iq«x  )  X(q,s  )+R(q)exp[-2iq(N-l)a  |cos  a  |]exp(iq^*x^) 


x  A(q,3)  ,      for  Zp  >  0  and  cos  a  <  0   .   III-33 


U5 


In  obtaining  Eqs .  III-33,  we  have  taken  the  phase  of  the  incident 
wave  to  be  unity  at  z  =  0. 

3.   Calculation  of  the  B^^(x. ) 

We  now  turn  our  attention  to  the  calculation  of  the 
-> 

B^  (xp)  which  must  be  used  in  the  evaluation  of  Eq.  11-30.   If 
Eqs.  III-33  are  substituted  in  Eq.  11-32,  one  finds  that  Eq.  11-26 
is  satisfied  provided  the  normalization  constant  A(q,s)  is  given  by 


A(q,s)  =  -^   .  111-3!+ 

>/nM 


In  evaluating  Eq.  111-3^,  we  have  noted  that  the  sum  in  Eq.  11-26 
may  be  performed  only  over  the  pure  regions  of  the  crystal.   This 
is  allowable  since  this  sum  produces  a  tenn  of  the  order  of  the 
volume  of  the  crystal.   In  the  limit  then  as  the  crystal  becomes 
infinite,  the  contribution  to  the  sum  for  points  x,  lying  within  the 
sample  becomes  negligible. 

We  therefore  have  from  Eq.  11-32  and  111-3^+  that  for  x^^- 
in  the  pure  regions  of  the  crystal 


B^^(x  )  =   ^  C^^(x^)     ,  III-35 

/n 


where  n  is  the  total  number  of  atoms  in  the  crystal.   For  future 
calculations  as  well  as  physical  interpretation  of  the  results,  it 
is  most  convenient  to  express  the  polarization  vectors  of  the  re- 
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fleeted  waves  as  a  function  of  the  propagation  vector  q  .   To  do  so 
we  adopt  the  following  convention:   One  transverse  polarization  vector, 
denoted  by  X(q,l),  will  be  directed  along  the  positive  x  axis.   (See 
fig.  1.)   Longitudinal  polarization  vectors,  denoted  by  A(q,3)  will  be 
directed  along  the  direction  of  propagation  of  the  wave.   The  second 
transverse  polarization  vector,  denoted  by  X(q,2),  will  lie  in  the 
plane  of  incidence  and  its  direction  defined  such  that 

t(q,2)  X   A(q,3)  -  A(q,l)   .  III-36 

Using  this  convention,  it  is  easy  to  show  that 

A(q,l)  =  X(q^,l) 

A(q,2)  =  -cos  2a  X(q^,2)  +  sin  2a  A"(q  3) 

t(q,3)  =   -cos  2a  A(q^,3)  -  sin  2a  A(q^,2)  .  III-37 

Use  of  Eqs.  III-37,  III-35  and  III-33  yields  the  following  rather 

->■ 

complicated  expressions  for  B^^Cx,): 
For  Zj  cos  a  >  0, 


B^^(Xj^)  =  -—  T(q)  exp[-iq(N-l)a|cos  a]] 
/n" 


exp(iq'x^)  A(q,s)   .  III-38a 


hi 


For  Zp  <  0,  cos  a  >  0, 


v<ll 


B   (x^)  =  [exp(ij»x^)  t(J,l)  +  R(J)exp(iq^'Xj^)  t(J^,l)] 

vn 


.12 


B   (x^)  =  [exp(iq  'x^)  A(q,2)  -  R(q)cos  2a  exp(ij^«x^)  A(J^,2) 


+  R(q)  sin  2a  exp(iq^»x^)  t(q^,3)] 


^q3_^ 

B   (x^)  =  [exp(iq'Xj^)  X(q,3)  -  R(q)cos  2a  exp(iq^»x^)  X(q^,3) 

vn 


-   R(q)  sin  2a  exp(iq^'Xj^)  A(J^,2)] 


III-38b 


For  z.  >  0,  cos  a  <  0, 


.ql 


B   (xj)  =  [exp(iq»x  )  A(q,l)  +  R(q)exp[-2iq(N-l)a|cos  a|] 


•n 


x  exp(iq^«x^)  X(q^,l)J 


.12 


B   (x^)  =  [exp(iq»x^)  A(q,2)  -  R(q)cos  2a  exp[-2iq(N-l)a|cos  a|] 


^ 


X  exp(iq^«x^)  X(q^,2)  +  R(q)sin  2a  exp[-2iq(lI-l)a|cos  a|] 
X  expdq^'x^)  A(q^,3)J 


.qs 


B   (xp)  =  [exp(iq»x.)  A(q,3)  -  R(q)cos  2a  exp[-2iq(N-l)a|cos  a|] 


v^ 
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X  exp(iq  'x^)  X(q  ,3)  -  R(q)sin  2a  exp[-2iq(N-l)a |cos  a|] 


->■  -)-  >  ->  ,-> 


X  exp(iq^v<j^)  A(q^,2)]  III-38c 


h.      Di 


scussion 


Eqs .  III-38  represent  the  central  result  of  this  chapter. 
In  view  of  the  vast  amount  of  mathematics  necessary  in  their  evalua- 
tion, as  well  as  the  fact  that  they  are  of  interest  quite  apart  from 
their  relevance  to  the  problem  at  hand,  we  interrupt  our  calculation 
of  the  heat  current  to  summarize  the  calculations  of  this  chapter  and 
to  discuss  the  physical  significance  of  the  results  obtained.   V/e  re- 
call  that  we  noted  that  the  quantities  B^  (x^, )  were  simply  related  to 
the  solutions  of  the  equations  of  motion  of  the  atoms  of  the  crystal 
which  could  be  written  as  the  sum  of  a  plane  wave  incident  upon  the 
sample  and  a  wave  scattered  by  the  sample.   We  began  by  considering 
a  single  plane  of  atoms  containing  isotopic  impurities  within  the 
crystal  and  found  that  such  a  plane  produced  a  reflected  plane  wave 
and  a  transmitted  plane  wave.   V7e  then  considered  all  possible  multi- 
ple reflections  of  these  waves  as  they  progressed  through  the  sample 
of  thickness  L   and  determined  the  total  wave  transmitted  by  and 

reflected  from  the  sample.   This  wave  was  then  used  to  determine  the 

-> 

quantities  B^^(x,)  given  by  Bqs .  III-33. 

We  begin  our  discussion  of  Eqs.  III-38  for  the  case  in 
which  2^  cos  a  >  0.   We  recall  that  a  is  the  angle  made  by  the  prop- 
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-> 
agatlon  vector  q  of  the  incident  wave  with  the  z  axis  and  that 

->• 

B^  (x»)  was  initially  written  as  the  sum  of  this  incident  wave  and 
the  wave  scattered  by  the  sample.   Thus  if  the  incident  wave  is 
incident  from  the  left  (cos  a  >  O)  and  if  one  observes  the  effect  at 
a  point  X.  for  to  right  of  the  sample  (z„  >  O),  one  finds  that  the 
scattered  wave  has  the  same  direction  as  the  incident  wave  and  com- 
bines with  it  to  produce  a  transmitted  wave.   The  polarization  of 
the  transmitted  wave  is  the  same  as  that  of  the  incident  wave  and 
consequently  produces  a  disturbance  in  the  same  direction  at  the 
point  Xp  as  would  the  incident  wave  in  the  absence  of  any  scattering. 
The  amplitude  of  the  transmitted  wave,  however,  has  been  reduced 
from  that  of  the  incident  wave  by  the  factor  given  in  Eq.  III-38a. 
This  reduction  in  amplitude  results  from  the  fact  that  part  of  the 
incident  wave  has  been  reflected  from  the  sample.   For  the  case  in 
which  both  cos  a  and  z-  are  less  than  zero,  exactly  analogous  re- 
marks apply . 

For  the  case  in  which  z.  cos  a  <  0,  the  situation  is 
slightly  more  complicated.   Eqs .  III-38b  predict  that  if  one  examines 
a  point  Xp  on  the  sane  side  from  which  the  incoming  wave  is  incident, 
one  observes  reflected  plane  waves.   These  waves  are  reflected  such 
that  the  angle  of  incidence  equals  the  angle  of  reflection.   If  the 
incident  wave  is  polarized  transversely  perpendicular  to  the  plane 
of  incidence  ('s  =  l),  there  is  only  a  single  reflected  wave  polariz- 
ed perpendicular  to  the  plane  of  incidence.   If,  however,  the  in- 
cident wave  is  polarized  longitudinally  (s  =  3)  or  transversely  in 
the  plane  of  incidence  (s  =  2),  the  reflected  wave  contains  compc- 
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nents  which  are  both  longitudinal  and  transverse  in  the  plane  of 
incidence.   These  remarks  are  illustrated  in  fig.  h   whore  the  vari- 
ous incident,  reflected  and  transmitted  waves  are  drawn  for  each 
type  of  incident  polarization. 

Clearly  the  problem  is  very  similaLr  to  the  reflection  and 
transmission  of  electromagnetic  waves  through  a  dielectric  interface. 
Tlie  only  significant  difference  is  that  lattice  waves,  unlike  electro- 
magnetic waves,  are  capable  of  propagating  longitudinal  modes.   It 
is  for  this  reason  that  one  finds  longitudinal  incident  waves  pro- 
ducing reflected  transverse  waves  as  well  as  reflected  longitudinal 
waves  and  vice  versa.   One  may  note  that  the  two  reflected  waves,  as 
does  the  transmitted  wave,  always  produce  a  net  disturbance  at  the 
point  X.  which  is  in  the  same  direction  as  a  disturbance  caused  by 
the  incident  wave.   Thus  if  the  incident  wave  is  transverse  perpen- 
dicular to  the  plane  of  incidence,  there  is  only  a  single  reflected 
wave  also  polarized  in  that  direction.   Similarly  incident  waves 
polarized  in  the  plane  of  incidence  do  not  produce  reflected  waves 
having  a  component  of  their  polarization  vector  perpendicular  to  the 
plane  of  incidence.   The  fact  that  a  wave  reflected  from  the  left- 
hand  side  of  the  sample  differs  by  a  phase  factor  from  a  wave  re- 
flected from  the  right-hand  side  results  from  the  fact  that  we  chose 
the  phase  of  the  incident  wave  to  be  unity  at  z  =  0.   This  differ- 
ence is  of  no  physical  importance. 

Before  concluding  this  chapter,  it  is  necessary  to  discuss 
the  approximation  of  replacing  the  sum  in  Eq.  III-2  by  the  integral 
of  Eq.  III-3.   To  do  so  we  noted  that  it  was  necessary  to  assume 
that  the  phonon  wavelength  was  much  longer  than  the  distance  between 
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(a) 


(b) 


%3 


q,3 


(c) 


Figure  4 

Directions  of  scattered  waves.   Waves  in  the  lower  left- 
hand  quadrant  represent  incident  waves,  those  in  the  upper  left-hand 
quadrant  reflected  waves  and  those  in  the  upper  right-hand  quadrant 
transmitted  waves.   The  wave  vector  and  polarization  index  of  each 
wave  have  been  indicated.   Each  of  the  three  possible  types  of 
polarization  of  the  incident  wave  is  represented  by  parts  (a),  (b) 
and  (c) .   Though  both  reflected  waves  are  reflected  at  angle  a, 
they  have  been  displaced  slightly  for  clarity  in  parts  (b)  and  (c). 
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impui'ities.   This  assumption  implies  that  the  phase  changes  negligibly 
in  Eq.  III-2  as  the  vector  x^^^  explores  the  impurity  plane  and  replace- 
ment of  the  sum  by  an  integral  should  be  a  reasonable  approximation. 
If  a  represents  the  fraction  of  atoms  in  the  plane  which  are  impurities, 
we  must  require  roughly  that 


X  »  2Traa~^/^   .  iii_39 


In  the  next  chapter  it  is  sho^ra  that  at  low  temperatures  this  condi- 
tion can  be  easily  satisfied  by  those  phonon  wavelengths  that  contri- 
bute significantly  to  the  heat  current. 


CHAPTER  IV 
HEAT  CURRENT 

In  this  chapter  the  heat  current  will  be  calculated  from 
the  solutions  to  the  scattering  problem  obtained  in  the  previous  chap- 
ter.  We  recall  that  those  solutions  were  valid  only  for  a  uniform 
distribution  of  impurities.   Therefore,  we  have  in  section  two  of  this 
chapter  obtained  an  approximate  expression  for  the  current  for  a  ran- 
dom distribution  believed  to  be  valid  at  somewhat  higher  temperatures. 
The  restrictions  placed  on  the  temperature  due  to  the  necessary  re- 
strictions on  the  phonon  wavelengths  in  the  solution  of  the  scattering 
problem  are  evaluated. 

1.   The  Uniform  Distribution 

Having  solved  the  scattering  problem  for  the  uniform  dis- 
tribution in  the  previous  chapter,  the  heat  current  may  be  evaluated 
by  Eq.  11-5^  in  which  we  found 


UkT  V    -^  ->,      s,s' 
s   q,q' 


-     lG_,(q,q')"' 


'=''  (e^-1)^ 


zss      zs  s        (e^-l)(e^  -l) 


IV-1 
6(aj(q,s)-w(q',s'))  AT^ 


where  G    ,(q,q')  is  given  by  Eq.  11-30.   The  calculation  of 


zss 


G    ,(q,q')  is  quite  tedious  and  adds  little  to  the  physical  under- 
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standing  of  the  problem.   We  have  therefore  performed  this  calcula- 
tion in  Appendix  B.   Before  indicating  the  results,  however,  it  is 
necessary  to  digress  briefly  in  order  to  indicate  one  important 
property  of  the  model  which  is  used  in  the  calculation  performed  in 
Appendix  B. 

Ordinarily  one  expects  to  have  to  use  Eq.  11-25  and  solve 
a  secular  equation  to  determine  the  normal  mode  frequency  spectrum. 

However,  this  necessity  can  be  avoided  in  the  present  case  by  use 

-30 
of  a  theorem  due  to  Lord  Rayleigh.    This  theorem  may  be  stated 

as  follows:   If  a  single  mass  M  contained  within  a  system  of  equal 
masses  connected  to  one  another  by  springs  obeying  Hooke's  Law  is 
reduced  (increased)  by  an  amount  61M,  all  frequencies  are  unchanged 
or  increased  (reduced)  but  never  by  an  amount  greater  than  the  dis- 
tance to  the  next  unperturbed  frequency.   The  model  at  hand  consists 
of  a  crystal  having  an  infinite  volume  and  thus  a  continuous  fre- 
quency spectrum.   Isotopic  impurities  are  contained  within  a  slab 
having  an  infinite  cross-sectional  area,  but  finite  thickness  L  . 
Consequently,  if  the  impurity  mass  M'  >  M,  the  introduction  of  any 
number  of  impurities  within  this  slab  cannot  cause  a  finite  shift 
in  the  normal  mode  frequencies.   The  same  is  not  true  for  M'  <  M, 
however,  as  some  of  the  normal  mode  frequencies  may  shift  to  values 
greater  than  those  for  the  pure  crystal.   These  may  be  shown  to  cor- 
respond to  complex  values  of  the  wave  vector  and  are  referred  to  as 

localized  modes.   At  low  temperatures  they  do  not  contribute  to  the 

19 
heat  current,    and  may  be  neglected.   V/e  can  therefore  use  Eq.  11-25 

with  M'  =  M  to  obtain  the  noi-mal  mode  frequencies.   We  obtain 
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p  ,k 


IV- 2 


Returning  to  the  evaluation  of  Eq.  IV-1,  we  note  that  it 
is  shown  in  Appendix  B  that  G    (q,q)  is  an  antisymmetric  function 
of  q.   Thus  the  first  tenn  in  brackets  of  Eq.  IV-1  does  not  contri- 
bute to  the  current.   In  addition,  it  is  shown  that  G   , (q,q' )  is 

s  o 
-*-->■->■-> 
non-zero  only  for  q'  =  q  or  q'  =  q  where 


In  =  q^i  +  q,J  -  q^k 


^1  =  "^x^   -^  S^^ 


From  Eq.  B-19,  we  find 


I    iG..(^,J')r 


s  .s 


ss 


'''o  ^r,   ^r, 

— —   CO  (q)  v^(q) 


|T(q)r  5-r,  +  |T(q)|2  lR(-)|2  5 


qq 


q^q 


IV- 3 


where  it  has  been  noted  that 

u(q,s)  =  (jL)(q) 

for  the  case  in  which  all  three  modes  of  polarization  propagate  at 
the  same  group  velocity.   In  Eq.  IV-3 


v(q)  =  v(q,s)  =  V(jj(q) 


IV- J+ 
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is  the  group  velocity  of  phonons  having  wave  vector  q,  L  is  the 
total  length  of  the  crystal  along  the  z  axis  (to  be  made  infinite 
in  the  final  foi-nulae),  and  T(q)  and  R(q)  are  to  be  obtained  from 
Eqs.  Ill- 31  and  III-32. 

As  noted  previously,  we  must  now  take  the  limit  as  the 

crystal  becomes  infinite  and  the  frequency  spectrum  becomes  con- 

->■     -^ 
tinuous.   Consequently,  sums  over  q  and  q'  appearing  in  Eq.  IV-1 

must  be  converted  to  integrals  using  the  prescription 


I      *     — ^  /  iq  IV-5 


3  "*■ 

where  V/(27t)   represents  the  density  of  states  in  q  space  and 

V  =  L  L  L  is  the  total  (infinite)  volume  of  the  crystal.   In  addi- 
X  y  z 

tion,  Kronecker  deltas  must  be  replaced  by  Dirac  delta  functions 
using  the  relationship 


5    •  ->   -^  5(q  -q')  IV-6 

q  q      L     ^x  ^x 

XX         X 


etc. 


becomes 


Using  Eqs.  IV-3,  IV-5,  and  IV-6,  we  find  that  Eq.  IV-1 


2 


T(q)r  6(q-q') 

IV-T 


+  |T(q)|^|R(q)|^  ^^V^J  '^^V'V^  '^^'lz''^J 


AT   . 


{e^-lf        ^ 
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In  obtaining  Eq.  IV-7,  it  has  been  noted  that  from  the  properties  of 
delta  functions 


5(q^-qJ  ^\-%^    6(co(q)-co(q'))  =    ^^   6(q-q') 

|v^(a)| 


for  q'  =  q 
and  IV- 8 


SCq^-q^)  6(q  -q  )  6  (co(q)-to(q' ) )   = 


|77^   5(q^-q;)  5(q^-q;)  SCq^nJ 


for  q'  =  q. 


The  integral  over  q  in  Eq.  IV-7  is  to  be  done  over  the  first 
Brillouin  zone. 

From  Eqs.  III-31  and  III-32,  one  finds 


|T(J)|2 


1  H-  x2(q)  U2_^(0) 


and  IV-9 

^  „    x2(J)  U^   (0) 


1   +  x2(q)  u2_^(0) 
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where  x(q)  and  0  may  be  obtained  from  Eqs .  III-12  and  III- 30,  re- 
spectively.  If  one  substitutes  Eqs.  IV-9  into  Eq.  IV-7  and  performs 

-> 
the  integral  over  q' ,  one  finds 

2  •  E 

"^•^z^  =  ITT  J  ^^  "^  ^^^  1^7^^^'  ~? 2 

^    (27t)\t'^  ^     (e^-1)^ 


IV-10 


. ^rp 

2  2  2  2^  z 

4COS  a 


where  explicit  use  has  been  made  of  Eq.  III-12.   In  the  acoustic 
approximation,  we  have 


2,^.  2  2 

w  (qj  =  V  q 


and 


|v  (J)|  =  v|cos  a]    .  IV-11 

Making  these  substitutions  in  Eq.  IV-7  and  ',«-iting  the  integral  in 
terms  of  the  dimensionless  variable 


^    kT 


one  finds  after  considerable  algebraic  manipulation  that 
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7T  a 


T 


.„/T 


T    I  ^      ^     — f ?   ^^ 

I 'lJ  0  (e^-l)2 


IV-12 


/    iLiz ^t' 

0            e  a^(T/T  )^C^         .    2   ,,,      -1  _. 
1   +  D  sm      (Ncos        Q) 


z 


2  2 

y  1-0 


T     is   the  Debye   temperature   for  this  model  which  has  been  defined   as 


-B  =  IeT 


and  y  -  cos  a.   Also  in  obtaining  Eq.  IV-i2,  we  have  made  explicit 
use  of  Eq.  III-26.   We  may  note  that  in  terms  of  the  parameter  C  > 


ea(T/T  )? 
0  =  cos[2(T/T^)Cy]  +  ^^—  sinC2(T/T  )Cy]  .    IV-II+ 


In  the  limit  of  low  temperatures ,  we  may  replace  the  upper  bound  on 

r   r   2 
the  C  integral  in  Eq.  IV-12  by  infinity  since  the  factor  e  /(e  -l) 

causes  the  integral  to  become  negligibly  small  for  large  values  of 

C   For  the  same  reason,  we  may  approximate  sin[2(T/T  )^y]  by 

2(T/Tj^)^y  and  consequently  0  by  cosC2(T/T  )Cy]  in  the  term  1-9^ 

yielding 


1  -  0^  =^  sin^[2(T/T^)Cy]  =  l+(T/Tj^)Vy^   .       IV-15 
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We  may  not,  however,  make  this  approximation  in  the  term  sin(Ncos~  0) 
without  further  consideration  since  this  term  is  oscillatory  and  con- 
tains in  its  argument  the  factor  N  which  can  be  quite  large. 
These  approximations  yield  finally  for  Eq.  IV-12 


-6kv 
7T  a 


T 


3  -  U    e^ 


/C 


1 
X  / 


0    (e^-l)2 


y^  dy 


dC 


P  ? 
0   h       e  o        .    2 
y    + 


AT 


-1 


r —  sin  (Ncos~  0) 


IV-16 


Unfortunately,  this  integral  cannot  be  evaluated  analytically.   How- 
ever, some  special  cases  may  be  considered. 

First,  if  e  =  0  the  integral  can  be  easily  performed  yield- 
ing 


<J  >„  = 


-72kv  /  T  ' 


z  0    2  3    T^ 
7T  a      D 


C(^)  AT 


IV-IT 


where  C  is  the  Riemann  zeta  function  and  we  have  used  the  fact  that 


n   e 


/  r  — f ^  d^  =  nl  C(n)  . 

0    ie^-lf 


'    iv-18 


For  large  n,  ^(n)  =  1.   Eq.  IV-l8  represents  the  heat  current  of  the 
isotopically  pure  crystal  and  agrees  with  expressions  obtained  pre- 
viously by  Allen  and  Ford   and  by  Erdos . 
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Second,  in  the  case  of  extremely  low  temperatures,  one  may- 
show  that  if  N  T/T  «  1,  the  term 


sin^(lIcos~^  0)  -   Un^(T/T  )^C^y^   . 


IV-19 


If  we  use  Eq.  IV-19  in  Eq.  IV-l6,  we  find 


<J>  =  =^ 

2     2   3 
7T  a 


/  rp   \3 


D 


/c 


0    {e^-lf 


d? 


1 
X  / 


3  , 

y  dy 


AT 


0  y^  +  NVa^(T/T^)V   "^^ 
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This  integral  may  be  approximated  for  small  ea   by  relatively  straight- 
forward though  extremely  tedious  means  to  produce 


<J  >  =^<J  >^  <^l  +  N^e^a^(T/T^)^ 
z     z  0   I  ^  '  D 

IV- 21 


li+T  +  6o(iog(Nea  T/T  )  -  y) 


where  <J  >q' is  given  by  Eq.  IV-IJ  and  y  is  Euler's  constant 


Y  =  0.5772 


Finally,  we  consider  the  more  likely  case  where  N  T/T  »  1. 
For  that  case  Allen  has  developed  a  procedure  for  approximating 
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integrals  of  the  form  of  Eq.  IV-I6  employing  integration  in  the 
complex  plane.   This  procedure  will  not  be  discussed  but  the  results 
of  its  application  to  Eq.  IV-I6  reveal  that  for  T/T^   »   1, 


<J  > 
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<J  >„ 
z  0 


1  + 


2  2 
e  o 


€0 
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A  far  simpler  though  less  rigorous  approximation  than  that  develop- 
ed by  Allen  may  be  employed  for  the  case  where  eo  «   1.   For  this 

2     —1 
case,  the  tenn  sin  (Ncos   0)  will  oscillate  rapidly  if  N  T/T  >>  1 

and,  for  small  CO,    it  should  not  be  a  bad  approximation  to  replace 

2     -1 

sin  (llcos   0)  by  its  average  value  of  1/2.   If  this  is  done  the 

integrals  in  Eq.  IV-I6  become  elementary  and  the  result  is 


<J  >  =  <J  >^ 

z      z  0 


1  - 


I4/2 
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which  is  in  close  agreement  with  IV-22  for  small  eo . 

Discussion  of  Eqs .  IV-IT,  IV-21  and  IV-22  will  be  deferred 
until  the  next  chapter.   Before  concluding  this  section,  however,  it 
is  perhaps  worthwhile  to  discuss  some  of  the  approximations  necessary 
in  the  solution  of  the  scattering  problem.   It  may  be  recalled  that 
for  the  solutions  of  the  scattering  problem  for  the  uniform  distri- 
bution to  be  valid,  we  had  to  require  roughly  that 


qa 
.1/2 


«  1 
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vhere  o  represented  the  fraction  of  impurities  in  any  given  plane. 
In  terms  of  the  variable  E,    ,    this  condition  becomes 


2(T/Tj^)5 

wj—  «  1   .  IV-25 

a 


It  may  be  noted  that  the  function  ?  e^/(e  -l)^  is  sharply  peaked 
and  attains  its  maximum  value  at  about  ^  =   h,    falling  off  rapidly 
for  larger  and  smaller  values  of  ^  .   Thus  the  condition  of  Eq.  IV-25 
may  be  replaced  again  approximately  by 


8  T/T 


or 


T       0^/2 

TfT-  «  —5-   .  IV-26 

D 


In  a  typical  dielectric,  T  is  of  the  order  of  several  hundred 
degrees.   Hence  at  low  temperatures,  Eq.  IV-26  can  be  satisfied 
quite  easily  even  for  relatively  small  values  of  a  and  we  conclude 
that  use  of _ Eqs .  III-38  has  been  appropriate  for  the  calculation  of 
the  heat  current.   As  a  final  point,  it  was  noted  in  Chapter  III 
that  the  current  tended  to  zero  as  cos  a  -^  0.   This  may  be  seen  by 
noting  from  Eq.  IV-l6  that  the  integrand  becomes  proportional  to 
y  as  y  =  cos  a  tends  to  zero. 
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2.      Approximation  for  the  Random  Distribution 

An  approximate  value  for  the  heat  current  for  the  case  in 
which  the  impurities  are  randomly  distributed  throughout  the  sample 
is  obtained  in  this  section.  It  is  assumed  that  the  calculation  of 
the  quantities  G    , (q,q')  may  be  performed  using  the  Born  approxi- 
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mation  for  the  B   (x^).   It  is  further  asstimed  that  the  phonon  wave- 
lengths are  not  long  compared  to  the  inter- impurity  separation,  i.e. 


-,   ^  2Tra 


where  a'  represents  the  fraction  of  impurities  in  the  entire  sample. 

Thus  for  the  results  to  be  valid  at  relatively  low  temperatures,  it 

must  be  assumed  that  the  fraction  of  impurities  is  quite  small. 

Using  these  approximations  we  show  in  Appendix  C  that  the 
term 


I  1    .  ^..„(J.^)G   ,..(^',q') 


q^qi   ^>^  Ve  --l;(e  -1) 
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X  6(a3(q)  -  w(q»  )) 


which   appears   in  Eq.    II-5I+   does   not   contribute   to  <J  >    in  the   lowest 

z 

non-vanishing  order  of  e.   In  addition,  it  is  shown  that 


65 


12L 
s.s'  L 


2  ^     £  a  La'  CO  (q)  v  (q) 
V  (q) 2 5_ 


hirv' 


IV- 29 


to  lowest  non-vanishing  order.   If  we  make  this  substitution  in 

Eq.  11-5^,  take  the  limit  as  the  size  of  the  crystal  becomes  infinite, 

and  proceed  along  steps  identical  to  the  development  of  Eq.  IV-l6,  we 

find 


^T  ^  -   -3TTh^    r    .^     2,-^.        e' 
<J  >  =  r — ^  J  dq  u  (q)  — p— 
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If  use  is  made  of  Eq.  IV-13,  Eq.  IV-30  may  be  written 
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As  T/T  ->  »°  ,  this  integral  may  be  easily  performed  yielding 
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Before  concluding  this  chapter  and  turning  to  a  general 
discussion  of  the  results,  we  shall  again  determine  the  approximate 
temperatures  for  which  one  might  expect  Eq.  IV-32  to  be  valid.   We 
initially  assumed  in  the  calculation  that 


qa  >  a 


,1/3 
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which  in  terns  of  the  variable  C  becomes 


2(T/T^)C  >  a 


.1/3 


IV-3^ 


Again  the  function  C  e  /(e  -l)   is  sharply  peaked  attaining  its 
maximum  about  E,   -   Q,      Thus  Eq.  IV-3^  becomes 


iGiT/T^)      >     a'^/^ 


which  implies 


T  > 
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Examination  of  Eq.  IV-32  reveals  that  the  second  term  in  brackets  may 
become  quite  large  at  termperatures  satisfying  Eq.  IV-35  except  for 
very  small  values  of  N,  e,  and  o'.   As  the  second  term  in  brackets 
approaches  unity,  one  may  suspect  that  the  Born  approximation  is  no 
longer  valid.   Thus  the  restriction  on  the  values  of  q  coupled  with 
restrictions  brought  about  by  the  Born  approximation  leads  us  to  con- 
clude that  this  treatment  for  the  random  distribution  may  be  expected 
to  be  valid  only  for  rather  limited  values  of  the  parameters  involved. 
This  point  will  be  further  discussed  in  the  next  chapter. 


chapter  v 
suj#iaey  and  discussion 

The  purpose  of  this  concluding  chapter  will  be  to  summarize 
the  calculations  made  and  to  attempt  to  discuss  the  motivation  which 
has  led  to  these  calculations.   The  plausibility  of  the  results  will 
be  discussed  with  reference  to  the  results  of  the  one-dimensional 
model  for  which  calculations  can  be  made  in  a  more  rigorous  manner. 
For  reasons  to  be  discussed,  the  results  here  preclude  detailed  com- 
parison with  experiment  or  even  with  other  theoretical  calculations. 
Nevertheless,  where  applicable,  attempts  will  be  made -to  explain  why 
the  results  are  reasonable  from  a  physical  standpoint. 

The  primary  emphasis  of  this  investigation  has  been  the 
formal  exposition  of  a  problem  in  three  dimensions  equivalent  to  that 
in  one  dimension  initially  set  forth  by  Allen  and  Fordl^  Though  we 
ultimately  hope  to  be  able  to  explain  experimentally  observed  phe- 
nomena, the  intent  of  this  dissertation  has  not  been  to  obtain  re- 
sults which  could  be  fitted  to  the  experimental  data.   Indeed,  the 
calculations  have  not  been  restricted  to  those  systems  for  which 
experimental  results  are  available.  Rather  we  have  sought  to  intro- 
duce some  general  techniques  for  treating  the  problem  in  three  dimen- 
sions and  hopefully  to  provide  some  insight  into  the  naLure  of  energy 
transport  in  three-dimensional  solids.   Thus  a  rather  detailed  cal- 
culation of  the  current  is  given  for  a  ciystal  containing  isotopic 
impurities  uniformly  distributed  throughout.   Tlie  results  of  this 
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calculation  are  given  by  Eqs .  IV-17,  IV-21  and  IV-22.   The  mathemati- 
cal techniques  used  in  this  solution  will  be  applicable  to  more  gener- 
al distributions  of  impurities  as  well.   This  point  will  be  further 
discussed  later.   In  addition,  an  approximate  solution  is  given  for 
the  case  in  which  the  impurities  are  randomly  distributed  throughout 
the  sample.   The  result  is  represented  by  Eq.  IV-32. 

If  one  wishes  to  define  a  thermal  conductivity  to  be  as- 
sociated with  these  currents,  one  may  do  so  by  writing 


<j  >   =   _K  1^ 
Z  dz 


V-1 


where  K  is  the  thermal  conductivity  and  8T/9z  the  temperature  gradient 
which  is  assumed  to  be  adequately  given  by 


8T  ^  AT 
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Using  this  definition,  the  thermal  conductivity  of  the  isotopically 
pure  three-dimensional  crystal  (see  Eq,  IV-lj)  becomes 
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For  the  case  in  which  impurities  are  uniformly  distributed  through- 
out the  sample,  one  finds  from  Eqs.  IV-21  and  IV-22, 


K  =  K 


I'- 


NVa^(T/T^)^ 


IU7  +  60(log(NGa  T/Tj^)  -  y) 
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for  N  T/Tj^  «  1,  and 
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for  N  T/T^  »  1.   Finally  for  the  case  in  which  the  impurities  are 
randomly  distributed  throughout  the  sanple,  we  have  from  Eq_.  IV-32, 


^o[i- 
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It  must  be  strongly  emphasized,  however,  that  the  values  of 
the  thermal  conductivity  obtained  in  Eqs .  V-3  through  V-6  all  depend 
on  the  thickness  L^  of  the  sample.   Consequently,  they  do  not  repre- 
sent thermal  conductivity  in  the  usual  sense  and  are  analogous  to  the 
results  obtained  in  one  dimension  by  Allen  and  Ford"*"^  mentioned  pre- 
viously.  Since  experimentally  solids  are  known  to  obey  Eq.  I-l  where 
K  is  independent  of  the  thickness,  some  explanation  of  these  results 
is  necessary.   The  variance  of  the  results  with  Eq.  I-l  may  be  ex- 
plained by  noting  that  the  Kubo  theory  obtains  an  expression  for  the 
current  which  may  be  written  as 


J   =  -K  1^  . 
z      9z 


V-7 


where  the  dots  represent  higher  order  terms  in  powers  of  the  gradi- 
ent.  Eq.  11-11  represents  the  first  term  in  this  expansion  and 
higher  order  terms  are  assumed  to  be  negligible.   Furthermore,  the 
gradient  which  appears  in  Eq.  V-J  is  that  which  is  initially  imposed 
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on  the  sample.   It  is   assumed  in  the  Kubo  theory  that  this  gradient 
is  not  significantly  altered  as  the  steady  state  is  approached.   It 
has  been  pointed  out  that  this  assumption  is  probably  reasonable  in 
systems  for  which  there  is  a  strong  scattering  mechanism.   Otherwise, 
however,  as  has  been  pointed  out  by  Allen  and  Ford,    the  Kubo  theory 
will  not  produce  a  thermal  conductivity  in  the  usual  sense.   Never- 
theless, one  may  still  use  Eqs .  IV-17,  IV-21,  IV-22  and  IV-32  to 
discuss  the  steady  state  current  in  tei-ras  of  the  temperature  differ- 
ence between  the  isotopically  pure  regions  of  the  crystal.   V7e  may 
conclude  from  the  results  above  that  multiple  scattering  effects 
from  uniformly  distributed  impurities  do  not  limit  the  flow  of  energy 
sufficiently  to  produce  a  thermal  conductivity  independent  of  the 
thickness  of  the  sample.   Some  additional  discussion  of  this  point  is 
given  later.   We  now  turn  to  a  discussion  of  specific  results. 

For  the  case  e  =  0,  one  merely  has  the  current  of  the  iso- 
topically pure  crystal  given  by  Eq.  IV-17.   If  one  notes  that  the 
specific  heat  of  the  phonon  gas  is  given  by 
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Eq.  IV-I7  may  be  written  in  the  form 


-3 

z    "T     z 
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This  equation  may  be  obtained  from  the  results  of  Allen  and  Ford 

22 
and  of  Erdos   by  setting  the  velocities  of  propagation  of  the  vari- 
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ous  polarization  modes  equal.   Eq.  V-9  is,  of  course,  non-divergent 
and  its  meaning  is  rather  obvious.   It  simply  states  that  in  the  ab- 
sence of  any  scattering,  the  current  is  ultimately  limited  by  the 
amount  of  energy  available  for  transport. 

In  considering  Eq.  IV-21,  we  may  note  that  for  extremely 
small  values  of  W  T/T  ,  £,  and  a   the  dominant  term  in  that  expres- 
sion is  the  one  involving  the  logarithm.   Thus  for  purposes  of  dis- 
cussion, we  write 


<j  >  ~  <j  > 

z  z  0 


1  +  60N^e^a^(T/Tj^)^  log(Nea  T/T^) 
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Since  Nea  T/T  <<  1,  the  second  term  in  brackets  is  negative  and  the 
current  is  reduced  from  that  of  the  pure  crystal  as  expected.   It 
might  at  first  appear  surprising  to  find  that  the  correction  term 
varies  as  {l/T^)      log(T/T  )  since  it  may  be  observed  in  Chapter  III 

that  the  intensity  of  a  wave  scattered  by  a  single  impurity  varies 

1+ 
as  0)  .   This  is  just  Rayleigh  scattering  and  one  would  ordinarily 

7 
expect  it  to  produce  a  term  proportional  to  T  as  has  been  found  in 

the  phenomenological  theory.   The  reason  for  the  weaker  dependence 

on  T  is  that  although  the  amplitude  of  a  wave  scattered  by  a  single 

2 
impujcity  does  indeed  vary  as  OJ  ,  the  total  amplitude  of  the  plane 

wave  scattered  by  an  entire  plane  of  impurities  is  a  more  compli- 
cated function  of  u.   For  small  frequencies  and  angles  not  too  close 
to  —  ,  the  amplitude  is  proportional  to  w.   This  difference  results 
from  the  fact  that  the  scattering  has  been  treated  coherently  with 
very  definite  phase  relationships  between  individually  scattered 
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waves  accounted  for.   Clearly,  this  is  the  proper  treatment  of  a 
scattering  problem  for  which  inpurities  are  uniforinly  distributed 
and  for  which  the  phonon  wavelengths  are  long  compared  with  the 

distance  between  successive  impurities. 

2 

The  dependence  of  the  correction  term  on  N  log  N  rather 

than  N  may  be  explained  similarly.   One  may  show  that  at  temperatures 
satisfying  the  condition  N  T/T  <<  1,  only  those  phonons  having  wave- 
lengths greater  than  the  thickness  L  of  the  sample  contribute  signi- 
ficantly to  the  current.   Thus  the  waves  scattered  by  the  planes  of 
impurities  within  the  sample  are  very  nearly  in  phase.   For  sources 

emitting  sound  waves  in  phase,  it  is  a  well-known  fact  in  the  theory 

2 

of  sound  that  the  resultant  intensity  is  just  N  times  greater  than 

that  from  any  one  source.   The  result  here  is  perfectly  analogous. 
The  weaker  dependence  on  log  N  results  from  the  fact  that  the  scat- 
tering has  been  treated  non-isotropically  and  is  rather  strongly 
dependent  on  angle. 

It  is  well  to  point  out  at  this  point  that  Allen  et  al. 
have  obtained  a  result  in  one  dimension  for  the  case  of  impurities 
uniformly  distributed  throughout  the  sample  that  is  very  much  ana- 
logous to  Eq.  IV-l6.   For  that  case  one  may  also  show  that  the  cor- 

2 
rection  term  to  the  current  is  proportional  to  N  ,  where  N  is  the 

3 

number  of  impurities  in  the  sample  and  to  T  .   The  difference  in  the 

temperature  dependence  results  from  the  difference  in  the  density 
of  states  at  low  frequencies  in  one  and  three  dimensions. 

Finally,  it  should  be  mentioned  that  although  Eq.  IV- 21 
was  derived  with  the  assumption  that  the  impurities  were  uniformly 
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distributed  throaghout  the  sample,  the  result  for  this  special  case 
(N  T/T  <<  l)  m.i^f  well  be  approximately  valid  for  more  general  dis- 
tributions of  impurities  as  well.   No  rigorous  proof  of  this  asser- 
tion will  be  given  but  it  can  be  made  plausible  by  a  physical  argu- 
ment.  For  phoncn  wavelengths  of  the  order  of  magnitude  of  the 
sample,  one  would  not  expect  that  the  deviations  of  the  impurities 
from  their  uniformly  distributed  positions  would  be  strongly  felt. 
This  results  frosi  the  fact  that  the  waves  scattered  by  individual 
impurities  will  be  very  nearly  in  phase  regardless  of  how  the  im- 
purities are  distributed  throughout  the  sample.   An  analogous  result 
has  been  noted  in  the  one-dimensional  model  vrhere  it  is  found  that 
the  current  obeys  an  equation  similar  to  IV-21  even  for  a  perfectly 
random  distribution  at  temperatures  such  that  N  T/T  <<  1  where  H 
is  the  number  of  impurities. 


For  the  case  in  which  N  T/T  «  1,  it  was  found  in  Eq.  IV-22 


that 
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This  expression  accounts  for  all  the  effects  of  multiple  scattering 
and  no  assumption  has  been  made  regarding  the  magnitudes  of  e  and  a. 
We  note  immediately  that  the  result  is  independent  of  the  number  of 
planes  N  and  thus  the  therrflal  conductivity  depends  linearly  on  the 

thickness  of  the  sample.   An  analogous  result  has  been  found  in  one 

2h  23 

dimension  by  Allen  et  al.    and  by  Greer  and  Ruben.     It  is  also 

3 

interesting  to  note  that  the  current  depends  on  T  AT,  a  dependence 
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39 
which  is  indicative  of  Kapitza  behavior   in  three  dimensions.   Thus 

if  one  employs  a  sample  in  which  all  atoms  have  mass  M' (a  =  l), 

Eq,  V-11  could  possibly  be  used  to  predict  the  experimentally  observed 

drop  in  temperature  as  the  current  proceeds  from  one  dielectric  to 

another  at  very  low  temperatures.   At  present,  the  model  has  not  been 

sufficiently  refined  to  account  for  all  the  effects  which  contribute 

to  Kapitza  resistance.   It  is  nevertheless  interesting  to  compare 

1+0 
qualitatively  Eq.  V-11  with  a  result  obtained  by  Little   who  has 

explicitly  treated  the  problem  of  Kapitza  resistance  in  solids. 
Little  considers  two  dissimilar,  semi-infinite  solids  joined  at  the 
plane  z  =  0  and  obtains  a  result  for  the  current  across  the  inter- 
face which  may  be  written  in  the  present  notation  as 

288  TTkV  C(^)  T^AT' 
<J  >  =  ^  „  ■  ■ „  V-12 

h  V 

At  is  the  temperature  difference  between  the  two  media  and  F  is  a 
complicated  integral  which  involves  the  total  transmission  coeffici- 
ent of  the  interface.   The  integral  is  not  explicitly  evaluated  by 
Little  but  he  shows  that  its  maximum  value  is  one-half.   The  maximum 
occurs  whenever  there  is  a  perfect  match  of  densities  at  the  inter- 
face.  If  one  assiimes  in  the  present  case  that  the  temperature  with- 
in the  sample  is  given  by  T  =(T  +  T  )/2,  then  the  temperature  dif- 
ference At'  between  the  isotopically  pure  region  on  th-^  left  and  the 
sample  is  given  by  AT'  =  AT/2.   Eq.  V-11  may  then  be  written  in  the 
form 
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If  one  further  notes  that  the  maximum  value  of  the  quantity  in 
brackets  in  Eq..  V-13  is  unity,  whereas  the  maximum  value  of  the  quanti- 
ty r  obtained  by  Little  is  one-half,  Eqs .  V-12  and  V-13  are  identical 
in  form.   The  examination  of  the  above  model's  ability  to  explain 
Kapitza  resistance  is  an  obvious  problem  for  future  investigation. 
In  fig.  5,  '^'^^>/<J^>Q   has  been  plotted  as  a  function  of 

the  product  ea.   It  is  seen  that  as  e  ^  ~  ,  <J  >  -»-  0.   This  results 

z 

from  the  fact  that  as  e  ^  oo  the  scattering  cross-section  becomes 
infinite  and  no  current  is  able  to  traverse  the  sample. 

Finally  it  is  interesting  to  discuss  Eq.  IV-32  with  regard 
to  the  one-dimensional  model.   It  is  seen  that  the  current  is  reduced 
from  that  of  the  pure  crystal  by  a  term  proportional  to  N  and 
(T/Tj^)  .   This  T  dependence  is  the  same  temperature  dependence  for 

the  correction  term  that  has  been  obtained  in  the  phenoraenological 

22 
theory  and  by  Erdos.    The  linear  dependence  on  the  number  of 

planes  results  from  the  fact  that  the  scattering  has  been  treated 
completely  incoherently  and  is  again  analogous  to  the  familiar  result 
of  incoherent  scattering  of  sound  waves  in  the  theory  of  sound. 
There  one  finds  that  if  the  phases  of  the  waves  emitted  by  the  vari- 
ous sources  are  com.pletely  random  the  resultant  intensity  is  N  times 

the  intensity  from  any  one  source. 

2^  23 

Allen  et  al.    and  Greer  and  Ruben   have  recently  succeed- 
ed in  obtaining  a  much  better  approximation  to  the  current  for  a 
randomly  disordered  chain  in  one  dimension.   Their  approximation  ac- 
counts for  multiple  scattering  and  is  valid  at  lower  temperatures. 
In  the  Born  approximation  of  their  expression,  one  also  finds  that 
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Figure   5 


<J  >/<J  >„  versus  ea  for  the  sample  conbaining  uniformly 


distributed  impurities. 
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the  current  is  reduced  from  that  of  the  pure  chain  by  a  term  which 
is  linearly  proportional  to  the  number  of  impurities  and  which  de- 
pends on  the  temperature.   The  temperature  dependence  is  different, 
of  course,  due  to  the  difference  in  the  density  of  states.   The 
author  suspects  that  Eq.  IV-32  may  well  be  valid  at  temperatures 
lower  than  those  assumed  for  the  purposes  of  making  the  calculation. 

It  does  not  seem  feasible  at  this  point  to  attempt  any 
comparison  of  the  results  obtained  here  with  experimental  data  for 
a  number  of  reasons.   First,  the  case  of  uniformly  distributed  iso- 
topes, for  which  a  rather  complete  calculation  has  been  given,  is 
except  for  the  case  of  Kapitza  resistance  clearly  not  relevant  to 
experiment.   Second,  the  approximation  obtained  in  Chapter  IV  for 
the  random  distirubtion  of  impurities  can  be  expected  to  be  valid 
only  for  very  limited  values  of  the  parameters  involved.   t-!ultiple 
scattering  effects  have  been  completely  neglected.   Finally,  it 
should  be  noted  that  so  far  as  the  author  is  aware,  the  only  experi- 
mental data  available  are  for  those  solids  for  which  one  would  expect 
boundary  scattering  to  diminate  at  low  temperatures,  that  is  for 
long  rods  having  small  diameters.   In  fact,  most  experiments  con- 
ducted at  low  temperatures  have  been  undertaken  expressly  to  verify 
the  effects  of  boundary  scattering.   These  effects  have  not  been 
considered  in  this  model.   Though  boundary  scattering  is  unquestion- 
ably important  in  solids  of  certain  geometrical  shapes,  the  author 
does  not  believe  that  it  is  of  fundamental  importance  in  understand- 
ing energy  transport  or  even  of  practical  importance  in  a  wide  variety 
of  solids.   It  is  hoped  that  experimental  evidence  on  solids  for 
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which  boundary  scattering  does  not  dominate  at  low  temperatures  will 
be  available  in  the  future. 

It  is  perhaps  well  to  summarize  the  results  of  Eqs.  IV-21, 
IV-22,  and  IV-32  "by  noting  that  in  each  case  the  current  has  been 
reduced  from  that  of  the  pure  crystal.   It  is  evident  that  a  formal- 
ism has  been  developed  from  which  one  can  study  the  effects  of  includ- 
ing scattering  processes  in  three  dimensions  which  has  been  the  in- 
tent of  this  dissertation.   Ultimately,  one  hopes  to  obtain  a  com- 
plete understanding  of  how  the  Fourier  heat  law  comes  about  in  three 
dimensions.   The  formalism  initially  set  forth  by  Allen  and  Ford  in 
1968  for  treating  the  problem  in  one  dimension  has  recently  shown 
considerable  promise  for  yielding  such  an  understanding  for  the  one- 
dimensional  case.  The   author  believes  that  a  similar  understanding 
is  possible  in  three  dimensions.   However,  the  mathematical  techni- 
ques in  three  dimensions  are  sufficiently  more  complicated  that 
working  out  the  details  may  require  some  time.   V/e  conclude  by  in- 
dicating some  specific  problems  needed  to  be  dealt  with  in  the  three- 
dimensional  problem. 

Of  foremost  importance  is  the  necessity  for  being  able  to 
treat  the  effects  of  multiple  scattering  for  more  general  distri- 

butions  of  impurities.   In  this  regard,  it  may  be  noted  that  Allen 

2I+ 
et  al.   have  just  recently  developed  a  scheme  for  treating  this 

problem  in  one-dimension.   In  essence,  they  have  found  that  the  Born 
series  may  be  reordered  in  such  a  way  that  the  deviations  of  the  im- 
purity positions  from  the  perfectly  uniform  distribution  may  be  aver- 
aged separately  from  the  sums  over  the  uniformly  distributed  posi- 
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tions.   Their  results  clearly  indicate  that  the  manner  in  which  the 
impurities  are  distributed  throughout  the  sample  is  of  utmost  impor- 
tance in  the  determination  of  the  thermal  conductivity  even  at  low 
temperatures.   Unfortunately,  such  a  procedure  has  not  yet  been  pos- 
sible in  three  dimensions  due  to  the  complicated  form  of  the  three- 
dimensional  Green's  function  and  the  fact  that  the  waves  scattered 
by  individual  impurities  are  spherical  rather  than  plane  waves.   In 
fact,  it  is  not  clear  that  one  can  always  obtain  plane  wave  solutions 
to  the  scattering  problem  for  an  arbitrary  distribution  of  impurities 
in  three  dimensions. 

The  second  important  feature  that  has  thus  far  not  been 
included  in  either  the  one-dimensional  or  the  three-dimensional  model 
is  the  important  effect  of  anharmonic  terms  in  the  crystalline  poten- 
tial which  give  rise  to  phonon-phonon  scattering.   It  has  been  assum- 
ed that  at  very  low  temperatures,  the  effects  of  these  terms  are  neg- 
ligible.  However,  at  high  temperatures,  the  lattice  motion  is  an- 
harmonic and  in  any  extension  of  this  model  to  the  high  temperature 
limit,  phonon-phonon  interactions  would  have  to  be  included. 


APPENDIX  A 

THE  METHOD  OF  STATIONARY  PHASE 
In  this  appendix,  we  seek  an  approximation  to  integrals  of 


the  form 


I  =   /   /  g(x',y')  exp[iqf(x',y')]  dx'  dy'  ,  A-1 

_oo  — oo 

The  approach  will  be  to  use  the  method  of  stationary  phase.   This 
method  has  been  discussed  for  the  case  of  single  integrals  by  vari- 
ous authors  and  comprehensively  treated  for  the  case  of  double  inte- 
grals when  f(x',y')  is  real  by  Jones  and  Kline   and  by  Kline  and 
Kay,    Since  the  rigorous  theory  is  quite  complicated,  our  discus- 
sion will  be  confined  to  some  rather  intuitive  arguments. 

If  the  function  qf(x*,y')  is  much  greater  than  unity,  one 
would  expect  the  exponential  term  to  oscillate  rapidly  and  tend  to 
average  to  zero.   The  only  place  where  this  would  not  be  true  is  in 
the  neighborhood  of  the  points  where  f(x',y')  is  stationary,  i.e. 
points  where  f(x',y')  obtains  a  maximum  or  minimiim  value.   Consequent- 
ly, a  good  approximation  to  the  integral  should  be  obtained  by  ex- 
panding f(x',y')  in  a  Taylor   series,  retaining  only  the  quadratic 
term,  and  integrating  the  resulting  expression.   The  points  at  which 
f(x',y')  is  stationary  can  be  found  by  requiring  that 


M(|l^  .  0    and    ill4^  =  0  ■  A-2 

dx'  dy' 
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and  solving  for  the  required  values  of  x'  and  y'. 

Denoting  such  a  point  by  (x  ,y  )  and  expanding  f(x',y') 
about  this  point,  we  have 


f(x',y')  =  f(xQ,yQ)  +|aQ(x'-XQ)^  +  |3Q(y'-yo)' 


1       '     ' 

+  -  YQ(x'-XQ)(y»-yQ)  +  ... 


A-3 


where 


^  _  9  f(x',y') 

'  ax.  2 


R  =  3  f(x\y') 


^O'^O 


8y' 


^O'^O 


Y. 


9  f(x',y') 
8x'  9y' 


t   I 


A-U 


Substitution  of  Eq.  A-3  in  Eq.  A-1,  yields 


oo  oo        r 
It         II.-        I  »  ' 

I  =  exp[iqf(x  ,y  )]  g(x  y  )  /  /  exp  }iq/2la^{x'-x^) 

— oo  — oo        I 


+  BgCy'-yQ)^  +  Yo(x'-XQ)(y'-yQ)][  dx'  dy' 


A- 5 


t   I 


where  we  have  assumed  that  g(x',y*)  can  be  replaced  by  gCx^jy.)  in 
the  integral.   This  should  be  a  valid  approximation  provided  g(x',y') 
varies  slowly  with  respect  to  the  exponential  as  it  does  for  the 
problem  at  hand.   The  integral  appearing  in  Eq.  A-5  can  be  integrated 
by  use  of  integrals  which  have  been  tabulated.    Performing  the 
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integral,  one  finds 


I  = 


27Ti         %         ,     '         \  r.   „/  •    • 


^J\^ 


^     sCxqjYq)  exp[iqf(xQ,yQ)]  A-6 


y;I 


where 


Qq  =  +  1    for    aQ3o>  Yo      a^  >  0 


Qq  =  -  1    for    aQ3Q>  Yq      a^  <  0 


Oq  =  -  i    for    aQ6o<  Yq 


In  the  most  general  case,  there  will  be  more  than  one  point  at  which 
f(x',y')  is  stationary  and  thus  I  becomes 


^  =  ~rl  ^  R(x'.,y.)  expriqf(x.,y.)]  A-7 

'-A/|a.3.-Y^| 


q   J  /: ~        J  J        J  J 


where  the  (x.,y.)  represent  points  at  which  f(x',y')  is  stationary. 
J   J 

For  the  problem  at  hand,  we  have  from  Eq.  III-2 

p         p    2  1  /? 
f(x',y')  =  [(x^-x')  +  y^-y')   +  Zj^]    +  y'  sin  a 

where  we  have  set  Xn  =  x' ,   Even  for  values  of  the  wavelength 
X  =   27T/q  which  are  many  lattice  spaces  long,  the  function  qf(x',y') 
is  much  greater  than  unity  except  as  both  x'  and  Xp  approach  the 
origin.   For  such  points  the  above  approximation  may  not  be  valid. 
However,  we  shall  assume  that  such  points  cannot  substantially  affect 
the  answer  and  use  the  approximation  of  Eq.  A-7  generally.   This 
point  has  been  further  discussed  in  Chapter  III. 


APPEiroiX  B 
CALCULATION  OF  G    , (q.q' )  FOR  THE  UNIFORM  DISTRIBUTION 


The  purpose  of  this  appendix  will  be  to  calculate  the 
^        ivq.,q.')  which  must  be  used  in  conjunction  with  Eq.  IV-1  in  the 
calculation  of  the  heat  current.   From  the  factor  6  (a)(q,s  )-oo(q' ,s  '  ) ) 
which  appears  in  that  equation  it  is  clear  that  our  only  interest 
is  in  the  case  where  Lo(q,s)  =  (jo(q',s').   In  Eq.  11-26  we  found  that 


^   o(x.-x.)   -^,  ,  _    *^ 

G_.^.(q,q')  =  I        I    (z.-z.)   «^   ^   J   b^  ^  ('x.  )B  ^^(x.)     &-1 

i , j  a , 3        /m.m. 
1  J 


•.ss.'J-')=..r  u.^...j:s^g^si'^\-iyfa,) 


and  have  previously  noted  in  Chapter  II  that  for  aj(q,s)  =  w(q',s') 
if  one  performs  the  sum  over  j  in  Eq.  B-1  the  resulting  expression 
is  independent  of   z.  .   Therefore,  if  we  fui'ther  assume  that 
*  o(x.-x.)  is  negligibly  small  when  x.-x.  is  larger  than  a  few 


n^'sV;r  ^  „.^  ^s, 


lattice  spaces,  we  may  use  the  values  for  B^        (x. )  and  B^  (x  )  in 

either  of  the  two  pure  regions  of  the  crystal  to  evaluate  Eq.  B-1. 

\Ie   begin  by  considering  the  case  in  which  both  cos  a  and 

cos  a'  are  greater  than  zero,  where  a  and  a'  are  the  angles  made  by 

->-     -> 

q  and  q'  with  the  z  a:<is .   For  this  case,  it  is  most  convenient  to 

choose  the  solutions  for  the  B^^(x.)  and  B^  ^  (x. )  far  to  the  right 

of  the  scattering  region,  i.e.  for  z.  and  z.  >  0.   Thus  from  Eq. 

J      1 

III-38a,  we  have 

^zss.^^'^')=i  X.  I  (--j^^g^v^j)^^^')^*^^) 

8U 
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X  exp[iq(N-l)a( |cos  a|  -|cos  a*|)]X  (q*,s') 


X   X  (q,s)  exp(iq'*x  )  exp(-iq«x.)  , 


By  making  a  change  of  variable  and  letting  x  =  x.-x.,  we  have 


B-2 


G_,(q,q')  =^     I         U^    ^,«(xj  T(q')  T  (q) 
i  ,m  a, 3 


X  exp[iq(N-l)a( |cos  a|  -  | cos  a' | )]  X  (q',s')  X„(q,s) 


X  expLi(q'-q) 'x.  )  J  expdq'x  )  .  B-3 


The  Slim  over  the  x  and  y  components  of  the  vector  x  can  be  done 

i 

immediately  by  using  the  fact  that 


y   exp[i(q-q' )*x. ]  =  n  n  5  •    5  •  B-h 

x.,y.  ^     ^^   ^xSc  W 

11 


where  n  and  n  represent  the  total  number  of  particles  along  the  x 
X      y   ^ 

and  y  axes.   Furthermore,  since  q  =  q'  and  both  cos  a  and  cos  a'  are 
greater  than  zero,  we  must  have  q  =  q  .   Thus  we  can  write 


y  exp[i(q-q')'x. ]  =  n^  6  '    6  '    5  '  B-5 

L        i-L  \H  H  /   ^j     0   q  q   q_  q_   q  q 

1  ^x^x   ^^   ^z^z 


where  n_  is  the  total  number  of  atoms  in  the  sample  of  length  L  , 


co2 
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For  further  simplification,  we  write  Eq.  IV-2  in  the  form 


^^'"^  ^ss'  =  Fa     ^,  ^a^^'^^   ^  ^aB^^m^  exp(iq.x  )  X  (q,s')   B-6 
a, 3         m 


where  we  have  used  the  fact  that 


I     A^(q,s)  X^(q,s')  =  6^^,   .  &.? 

a 


By  taking  the  gradient  of  expression  B-6  with  respect  to  q  one  can 
show  that  for  s  =  s ' 


^^  u(J,s)  =  v(J,s)  =  i:; I     X  (J,s) 

^  2Maj(q,s)   a,B 

B-8 
m 


where  v(q,s;  is  the  group  velocity  of  phonons  having  wave  vector  q 
and  polarization  s.   In  addition.  Hardy  has  shown  for  s  5^  s',  but 
w^q,s)  =  oo(q,s')  which  is  true  in  the  present  case  that 


^,   ^a^^'")  ^  W\^  \   ^^P(i<l-^„i)  X3(J,s')  =  0  .     B-9 
OtjP  m 


Use  of  Eqs.  B-5,  B-8  and  B-9  in  Eq.  B-3  yields 


— 2iL 

G^ss'^^'^')  =  ^(^(q,s)  v^(q,s)  |T(q)  |  ^  6^,  6^^ ,      B-10 

z 
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where  L  is  the  total  length  of  the  crystal  along  the  z  axis  (to  be 
made  infinite  in  the  final  formulae).   Eq.  B-10  was  derived  with  the 
assumption  that  cos  a   and  cos  a'  were  greater  than  zero.   However, 
if  both  cos  a  and  cos  a'  are  less  than  zero,  the  same  results  are 
obtained.   We  conclude  therefore  that 


2iL 

z 

B-11 
for  cos  a  cos  a'  >  0  . 


The  case  in  which  cos  a  >  0,  cos  a'<  0  is  slightly  more 

complicated.   For  the  first  of  these  cases,  we  choose  the  solutions 

-»■  -> 

for  B^  ^  (x.  )  and  b'^^(x.)  far  to  the  left  of  the  scattering  region, 

i.e.  z. ,  z  <  0.   We  can  write  from  Eq.  III-38b  the  general  expres- 
1   J 

sion 


*qs 


1 


B    (x.)  =  —  [exp(-iq-x.)X  (q,s)  +  f, (q,s)  exp(-iq  'x  ) 
pj/r         Jp        -L  -J-J 


y^ 


X  Ag(q^,s)  +  f^Cq.s)  exp(-iq^«x  )  Xg(q^,3)  5^^ 


*  ,-y 


+      f2(q,s)  exp(-iq^'x  )  X^(q^,2)  6^^] 


B-12 


where  the  f's  are  to  be  obtained  from  III-38b.   Similarly,  we  have 
from  Eq.  III-38a  for  this  case 


->• 
q's* 

B     (x. )  =  —  [T(q')exp[-iq(N-l)a|cos  a' |  ] 
ot     1      / — 
vn 
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X  exp(iq''x.)  X  (q',s')  ]    .  B-13 


By  proceeding  in  a  manner  identical  to  the  evaluation  of  B-11,  one 
finds 


2iL 
G    ,(q,q')  =-7-^  u(q,s)v^(q,s)T(q)  exp[-iq(N-l)a |cos  a]] 

ZS  S  J-i  z 

z 


for  cos  a  >  0,  cos  a'  <  0  B-lU 

where  we  have  used  the  fact  that 


^z^^l^  =  -''z^^^ 


and  that 


T(q)  =  T(q^) 


as  can  be  seen  from  Eq.  III-31. 

Finally,  for  the  case  in  which  cos  a  <  0,  cos  a'  >  0,  we 
again  choose  the  solutions  for  to  the  left  of  the  scattering  region 
and  proceed  as  before  to  obtain 

— 2iL 
G^^^,(q,q')  =  -;: a)(q,s  )v  (q,s  )T(q)  exp[-iq(N-l)a  |cos  a|] 

ZS  S  Li  Z 

z 
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X  [f^Cq'.sMS^,,  -   f2(^'>^')^'2^s3  "  ^3^^''^')^'3^s2^  ^J' 


q-q 


for  cos  a  <  0,  cos  a'  >  0  .  B-15 


We  see  immediately  from  Eqs.  B-11,  B-lU,  and  B-15  that 

.->■  -►  .  ->•->-♦•->■ 

G    (q.,q.' )  is  non-zero  only  for  q'  =  q  or  q'  =  q,  .   Furthermore, 
zss'  1 

we  see  from  Eq.  B-11,  that  G   (q.,q)  is  an  antisymmetric  function  of 
q  since 


V  (q)  =  -  V  (-q)   .  B-16 


Consequently,  the  first  term  in  hrackets  in  Eq.  IV-1  does  not  con- 
tribute and  we  need  only  evaluate 

I  lG^^g(J.5)|^  and   I         |G^33,(q,qi)|^   . 
s  s,s' 

From  Eq.  B-11,  we  find 

2 

I         |G^33,(q,q)|'=  ^a32(^)v'(q)  |t(J)|^         &-1T 
s,s'  L 


where  we  have  noted  to(q,s)  =  u(q)  and  v  (q,s)  =  v  (q)  for  the  case 


90 


where  v  =  v^  =  v.   In  a  similar  fashion,  by  obtaining  the  values  of 
the  f's  in  Eqs.  B-li+  and  B-15  from  Eq.  III-39b,  one  finds 


2 

|2    "^^^0   2,-^s   2,->^  i„,-^vi2  i„,-^m2 


I    l^zss'^'l'^l)!   =   -^u;^(q)v  (q)  |T(q)|^  |R(q)|' 


s,s'  L 

*  z 


or  in  general 


B-18 


2 

12L         ? 

I    |G^33,(q,q')|'=  -^O^^C^)  V  (q)   [lT(^)|^ 


X  6JJ,  +  |T(ci)l^  |R(<l)r  5-.-]    .       B-19 


^1 


APPENDIX  C 


APPR0Xir4ATI0N  OF  G    ,(q,q')  FOR  THE  RAITDOM  DISTRIBUTION 


We  sheill  obtain  an  approximation  for  the  quantity 

G    i(Q.»q.')  in  this  appendix  for  the  random  distribution.   The 

calculation  will  involve  the  use  of  the  Born  approximation  for  the 

solutions  of  the  B^  (xp)  and  a  random  phase  approximation.   From 
Eq.  11-26, 


^zss'^^'^') 


1    I  (--l) 

i,j  oc,3     ^ 


a3  i  '  .] 

•^        v-mTm". 
1  J 


(x.-x . )   q's ' 


a 


*qs 

U)  B    (x  ) 
1    tS     J 


C-1 


where  the  B   (x,  }  are  given  in  the  Born  approximation  by 


^    /n 


2  3    exp[iqx(in^, )] 


expdq.x.  )X(q.s)  -  — ^  I  -jjj^ 


Uttv   r 


T 


exp(iq'Xn^)  X(q,s) 


C-2 


If   Eq.    C-2   is   substituted   in  Eq.    C-1,   one   finds 


zss 


.(q,q')=i      I'       I    (z.-z    )    ""gj^i^i      >^(q.,s')    A    (q,s) 
i,j   a, 3  /m . m . 


exp[-iqx(jn^)] 


2  3 

expUq' 'x.  jexp(-iq*x. )   -  5—  exp(iq'«x.)   p r 


J-        l.r  2 


+TTV 


->.  ^  c-,,2^3  .    .        exp[iqx(in    )] 

X      exp(-iq'Xn^)   -  —  exp(-iq«x    )   ,         . exp(iq'-Xn^) 

I47TV  ^  ^      r' 
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C-3 

e^oo^a^  r   exp[iqx(inj.)]exp[-iqx(jn^,)]     ,  h-  -v   ,      ^^   '  1 
T~r  )   r- ^ r~' ■; exp(iq'*Xri  )exp(-iq*Xn  ,) 


For  a  random  distribution  of  impurities  and  for  cases  in  which 
qAxn  >  1,  one  would  expect  the  second  and  third  terms  in  brackets 
in  Eq.  C-3  to  contribute  negligibly  due  to  the  presence  of  random 
phase  factors.  Similar  statements  would  apply  to  the  fourth  term 
except  for  cases  in  which  r  =  r'  and  q  =  q*.   Thus  Eq.  C-3,  becomes 


i,j  a, 6     -^   /m^m 


exp(iq' 'x. ) 


C-U 


e^co^a^   r  exp[iqx(inj.)]exp[-iqx(jn^)] 


x  exp(-iq«x.)  +  A-  I   n \ r \ 6"^, 


The  first  term  in  brackets  can  be  shown  to  produce  just  the  value  of 

0 

G   __ ,  for  the  isotopically  pure  crystal.   Denoting  this  by  G    ,  and 
^S  ^  z*s  s 

obtaining  its  value  from  the  results  of  Appendix  B,  one  finds 


0    ^  ^     -2iLQ   ^ 

z  ^^ 


Denoting  the  contribution  from  the  second  term  in  brackets  in  Eq.  C-U 
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$  „(x.-x.)    _      ^     ^2  1+6 

^zss'^^'^  )  =  n  J.      ^«^V^j)   ,, X^(q,s)A  (q,s')  —j^ 

1  j 


X  )^  exp[iqx(in^)]  exp[-iqx( jn^.)]  g_^  j,_g 

r  r: c j-. r qq' 

x(in^)  x(jn^; 


For  the  case  in  which  cos  a  >  0,  we  choose  the  solutions  far  to  the 

left  of  the  sample  to  evaluate  Eq.  C-6.   Then,  noting  that  if 

$  o(x.-x.)  is  negligibly  small  unless  x.-x  is  less  that  a  few  lattice 
Ctp   1   J  1   J 

spaces,  we  can  write 


exp[iqx(in  )]  exp[-iqx( jn^)]  -  exp(iq'*x^-x  ) 


where 


»  =       ^    =    J   "r 


Then,  upon  setting  x  =  x.-x.,  Eq.  C-6  becomes 
»  -^        ^  m    1  J » 


0^3,, (q.i-)  =  5^  r  l^^  ^^aj  x^(j,.)  Xg(j.s') 

i,m  a, 3 


9U 


X 


e 


2  h  6  exp(iq' -x  ) 


I  —, .  C-T 


IDTT  V   r   X  (.in  J 
r 


As  in  Appendix  B,  it  is  now  necessary  to  express  the  polari- 
zation vectors  as  a  function  of  the  wave  vector  q'.   If  we  adopt  the 
convention  of  Chapter  III,  clearly  these  relationships  are 

X(q,l)  =  X(q',l) 

X(q,2)  =  -  cos  (a+6)  t(J',2)  +  sin  (a+6)  t(J',3) 

t(J,3)  =  -  cos  (a+B)  A(J',3)  -  sin  (a+B)  X(5',2)  .   C-8 

By  using  these  relationships  and  Eqs.  B-8  and  B-9,  it  can  he  shown 
that  Eq.  C-T  becomes 


2  6        o      5^^N 

^     z-^-  ->,  V    -1   V   £  a  V  cos  B  CO  (q;  r.     r       n  n 

i,r  OTT  V     X  (m  ) 
'  r 


The  sum  over  i  in  Eq.  C-9  may  be  converted  to  an  integral  since 

X  (in  )  is  a  large  number  which  changes  negligibly  as  Xn  changes, 
r  ^ 

Thus  by  noting  that 


Z.-Zn 

o  1  "r 

cos  B  =  —r r 

X  in  ) 
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we  have 


cos  6     .  a'   r  ^i"^^: 


2,.   ^       3   ^   3/.   ^ 
X  (in  /      a     X  (in^) 

r  T 


/  -^ ^    dxn^  C-10 


where  a'  denotes  the  density  of  impurities  in  the  sample.  Since  the 

impiirity  plane  is  infinite  the  integral  appearing  in  Eq.  C-10  cannot 

possibly  depend  on  the  x  and  y  components  of  the  vector  x. .   Thus 

they  may  be  set  equal  to  zero  for  the  piirpose  of  performing  the 
integral  and  one  easily  finds 


-  J^^r-T— .   =  —3—    •  '^-^^ 

a         X  (inj.;       a 


This  result  may  be  substituted  into  Eq.  C-9  and  the  sum 
over  i  performed  to  produce  the  number  of  atoms  in  the  sample.   The 
final  result  becomes 

?  "?  2   c  -^ 

ie  a^L„a'co-^(q) 

G    ,(q,q')  =  T- ^-^t  ^   . 

z 

for  cos  ex  >  0   ,  C-12 


For  the  case  in  which  cos  a  <  0,  one  may  evaluate  G   i(q,q')  ^Y 

zss 

choosing  solutions  for  the  B^  (x,)  far  to  the  right  of  the  sample. 
In  this  case  one  finds 
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4ttv  L 


C-13 


Combining  Eq.  C-5  with  Eqs.  C-12  and  C-13,  one  finds 


-2iLQa)(q) 


V  (q)  - 


2  ?     U,->x 
eVLQO'o)  (q) 


for  cos  a  >  0 


C-lU 


and 


^zss'^^'^') 


-2iLQw(q) 


qq'   ss ' 


v^lq 


(q)  ^  e^a^gCo^  (q) 


8ttv- 


for  cos  a  <  0 


C-15 


These  values  may  now  be  used  to  evaluate  the  quantity  in 
brackets  in  Eq.  II- UT.  One  finds  that 


I        G_(q,q)  G^  .  .(J',J') 


zss       zs ' s 
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does  not  contribute  to  the  current  in  the  lowest  non-vanishing  order 
of  e.   By  using  Eqs.  C-li+  and  C-15  one  may  also  show  that 


I         |G,..,(q,q') 


2  ?  ^. 

12L^W    (q) 


s,s' 


zss 


■  2^         e^L  a'a)\q)  |v   (q) 
V    (q)    -   


i+TTV" 


qq» 


C-l6 


to  lowest  non-vanishing  order  in  e. 
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